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Multiscale  Analysis  of  Heterogeneous  Media  in  the 
Peridynamic  Formulation* 

Bacim  Alali^  Robert  Lipton^ 

October  28,  2008 


Abstract 

A  rigorous  multiscale  method  is  presented  for  modeling  the  dynamics  of  fiber- 
reinforced  composite  structures  using  the  peridynamic  formulation.  The  multiscale 
analysis  delivers  a  new  multiscale  numerical  method  that  captures  the  dynamics  at 
structural  length  scales  while  at  the  same  time  is  capable  of  resolving  the  dynamics  at 
the  length  scales  of  the  fiber  reinforcement.  The  new  numerical  method  is  able  extract 
this  information  at  a  cost  that  is  anticipated  to  be  far  less  than  the  direct  numerical 
simulation  of  structural  components  made  from  multiple  plys  containing  thousands  of 
fibers. 
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1  Introduction 


The  peridynamic  formulation  introduced  in  m  is  a  nonlocal  continuum  theory  for  de¬ 
formable  bodies  that  does  not  use  the  spatial  derivatives  of  the  displacement  held.  In¬ 
teractions  between  material  particles  are  characterized  by  a  pairwise  force  held  that  acts 
across  a  hnite  horizon,  see  Section  |1.1[  The  same  equations  of  motion  are  applicable  over 


the  entire  body  and  no  special  treatment  is  required  near  or  at  defects.  These  properties 
make  it  a  powerful  tool  to  model  problems  that  involve  cracks,  interfaces,  and  other  defects, 
see  [21  m  [m  [22l  [23l  |2l] .  This  work  focuses  on  the  multiscale  analysis  of  heterogeneous  media 
using  the  peridynamic  formulation.  The  objective  is  to  provide  numerical  methods  that  cap¬ 
ture  the  dynamics  inside  composites  at  both  the  structural  scale  and  the  microscopic  scale 
with  a  cost  far  below  that  of  direct  numerical  simulation. 

We  consider  particle  or  hber  reinforced  composites.  Here  the  characteristic  length  scale 
of  the  particle  or  hber  reinforced  geometry  is  assumed  to  be  very  small  relative  to  the  length 
scale  of  the  applied  loads.  The  length  scale  of  the  microstructure  is  denoted  by  e.  We 
study  three  peridynamic  models  of  hber-reinforced  materials.  In  the  hrst  model,  which  we 
call  “the  short-range  bond  model”,  the  peridynamic  horizon  is  of  the  same  length  scale  as 
that  of  the  microstructure  and  the  horizon  approaches  zero  as  £  goes  to  zero.  In  the  second 
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Figure  1:  Fiber-reinforced  composite. 


model,  a  long-range  ^-independent  pairwise  force  is  added  to  the  short-range  pairwise  force 
of  the  hrst  model.  Here  the  long-range  pairwise  force  depends  only  on  the  relative  position 
of  the  two  particles  and  the  associated  peridynamic  horizon  is  hxed  and  independent  of  the 
microstructure  length  scale  £.  We  will  refer  to  the  second  model  as  “the  short-range  and 
long-range  bond  model”.  In  the  third  model,  we  consider  a  long-range  pairwise  force  that 
fluctuates  with  the  microstructure.  The  peridynamic  horizon  in  this  model  is  hxed  and 
independent  of  e.  This  model  will  be  called  “the  huctuating  long-range  bond  model”.  In 
all  of  these  models,  the  peridynamic  initial  value  problem  is  a  partial  integro-differential 
equation  with  rapidly-oscillating  coefficients  supplemented  with  initial  conditions. 

For  the  hrst  two  models  the  concept  of  two-scale  convergence,  introduced  by  Nguetseng 
[I8]  and  Allaire  [Tj,  is  used  as  a  tool  to  identify  both  the  macroscopic  and  microscopic  dy¬ 
namics  inside  the  composite.  A  downscaling  method  obtained  through  the  use  of  Semigroup 
theory  provides  a  strong  approximation  for  capturing  the  mirco-level  huctuations  about  the 
macroscopic  displacement  held.  The  multiscale  approximation  obtained  for  the  hrst  two 
models  are  shown  to  be  good  approximations  to  the  actual  solution  in  the  norm  when 
the  microstructure  is  sufficiently  hne.  Explicit  error  estimates  are  provided  for  sufficiently 
regular  initial  and  loading  data  for  the  hrst  model.  This  multiscale  analysis  provides  the 
theoretical  framework  for  a  new  multiscale  numerical  method  for  computing  the  deformation 
of  hber-reinforced  composites  in  the  presence  of  residual  forces.  The  multiscale  numerical 
method  delivered  here  captures  the  dynamics  at  structural  length  scales  while  at  the  same 
time  is  capable  of  resolving  the  dynamics  at  the  length  scales  of  the  hber  reinforcement.  The 
new  numerical  method  is  able  extract  this  information  at  a  cost  that  is  anticipated  to  be  far 
less  than  the  direct  numerical  simulation  of  structural  components  made  from  multiple  plys 
containing  thousands  of  hbers. 

For  the  third  model,  the  Semigroup  theory  of  linear  operators  [la  HB]  is  utilized  to 
identify  both  the  macroscopic  and  microscopic  dynamics  of  the  composite.  These  are  used 
to  develop  an  approximation  to  the  actual  solution  that  is  shown  to  be  a  good  approximation 
to  the  actual  solution  in  the  norm  when  the  microstructure  is  sufficiently  hne.  Explicit 
error  estimates  for  the  approximation  are  provided  for  this  model.  Last,  the  corresponding 
multiscale  numerical  scheme  is  presented. 

This  report  is  organized  as  follows.  Section  1^  provides  an  overview  of  the  peridynamic 
formulation  of  continuum  mechanics.  In  Section  1^,  we  introduce  three  peridynamic  models 
of  hber-reinforced  composites.  The  results  for  the  hrst  two  models  are  discussed  and  derived 
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in  Sections  |2]j^  In  Section  we  present  a  multiscale  analysis  method  for  these  two  models. 
Section  provides  uniqueness  and  existence  results  for  the  linear  peridynamic  initial- value 
problem  ( 1.10 )-( 1.12).  In  Section]^  we  review  two-scale  convergence  and  then  use  it  to 
identify  the  two-scale  asymptotic  limit  of  ( 1.10 )-( 1.12).  In  Section!^  we  build  on  the  analysis 
provided  in  Section  to  justify  the  results  of  Section  Section^  is  devoted  to  the  third 
peridynamic  model  of  hber-reinforced  composites.  A  multiscale  analysis  method  is  presented 
and  justihed  for  this  model. 


1.1  The  Peridynamic  Formulation  of  Continuum  Mechanics 

In  the  peridynamic  theory,  the  time  evolution  of  the  displacement  vector  held  u,  in  a  het¬ 
erogeneous  medium,  is  given  by  the  partial  integro-differential  equation 


p{x)  d^u{x,t) 


f{u{x,t) 


u{x,  t),x  —  x,x)  dx  +  b{x,  t), 


{x,  t)  E  fl  X  (0,  T) 


(1.1) 


where  is  a  neighborhood  of  x,  p  is  the  mass  density,  6  is  a  prescribed  loading  force  density 
held,  and  is  a  bounded  set  in  Here  /  denotes  the  pairwise  force  held  whose  value  is 
the  force  vector  (per  unit  volume  squared)  that  the  particle  at  x  exerts  on  the  particle  at 
X.  For  a  homogeneous  medium  /  is  of  the  form  f{u{x,t)  —  u{x,t),x  —  x),  i.e.,  it  depends 
only  on  the  relative  position  of  the  two  particles.  We  will  often  refer  to  /  as  a  bond  force. 
Equation  (1.1)  is  supplemented  with  initial  conditions  for  u{x,  0)  and  dtu{x,  0).  For  the  sake 
of  simplicity,  we  assume  constant  mass  density  given  by  p{x)  =  1.  However,  the  removal  of 
this  hypothesis  presents  no  barrier  to  the  subsequent  analysis.  For  the  purposes  of  discussion 
it  will  be  convenient  to  set 

f  =  X  —  X, 


which  represents  the  relative  position  of  these  two  particles  in  the  reference  conhguration, 
and 

rj  =  u(x,  t)  —  u{x,  t), 

which  represents  their  relative  displacement  (see  Figure]^.  In  the  peridynamic  formulation, 
it  is  assumed  that  for  a  given  material  there  is  a  positive  number  6,  called  the  horizon,  such 
that 

=  0,  for  1^1  >  6. 


The  pairwise  force  held  /  is  required  to  satisfy  the  following  properties: 

f{-v,  x  +  0  =  -fiv,  x)  (1-2) 

which  assures  conservation  of  linear  momentum,  and 

(^  +  h)  X  =  0 

which  assures  conservation  of  angular  momentum. 

A  material  is  said  to  be  microelastic  if  the  pairwise  force  is  derivable  from  a  scalar 
micropotential  oj 


=  —{r],f,x). 
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Figure  2:  Deformation  of  a  bond  within  the  peridynamic  horizon. 


It  can  be  shown  that  for  a  microelastic  material  the  pairwise  force  is  of  the  form  (see  lai) 

where  if  is  a  real- valued  function.  Finally,  a  material  is  linear  if  the  associated  bond  force 
f{ri,^,x)  is  linear  in  p. 

In  this  treatment,  all  materials  will  be  taken  to  be  microelastic  and  linear. 

1.2  Three  Peridynamic  Models  of  Fiber-Reinforced  Materials 

To  £x  ideas,  we  consider  a  periodic  medium  of  unidirectional  fiber-reinforced  material.  Here 
the  pairwise  force  is  given  by  the  linearized  version  of  the  bond-stretch  model  proposed  in 

m 

Here  a  is  a  real-valued  function  satisfying  a{x,x)  =  a{x,x).  We  will  study  three  different 
peridynamic  models  for  this  composite.  These  models  are  distinct  in  the  way  the  coeffi¬ 
cient  a  and  the  neighborhood  set  are  defined.  We  start  by  providing  the  mathematical 
description  of  the  periodic  microgeometry. 

Let  Y  C  be  a  unit  cube  and  the  local  coordinates  inside  Y  are  denoted  by  y  with 
the  origin  at  the  center  of  the  unit  cube.  The  unit  cube  is  composed  of  a  hber  which  is 
surrounded  by  a  second  material  called  the  matrix  material,  see  Figure]^  Let  Xf  denote  the 
indicator  function  of  the  set  occupied  by  the  hber  material  and  Xm  denote  the  the  indicator 
function  of  the  set  occupied  by  the  matrix  material.  Here  Xf  is  given  by 

/  \  _  j  y  is  in  the  hber  phase, 

Xi[y)  “  I  0,  otherwise, 

and  Xm  is  given  by 

Xm{y)  =  l-Xi{y)- 

We  extend  the  functions  Xf  and  Xm  to  by  periodicity.  For  future  reference,  we  denote  by 
6{  and  6^  the  volume  fractions  of  the  hber  material  and  the  matrix  material,  respectively. 
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Here  6*f  =  JYXf{y)dy  and  6*m  =  1  —  Of-  Also,  we  let  n  denote  a  unit  vector  parallel  to  the 
fiber  direction. 

In  the  first  model,  the  short-range  pairwise  force  is  given  by 


/short  iVyy  2/) 


oi{y,y  +  Qh^^Vy,  Kr/I<^ 

Isyl 


0,  otherwise, 

where  y  &  Y ,  =  y  —  y,  rjy  =  u{y,  t)  —  u{y,  t),  and  a  is  given  by 

a{y,  y)  =  Q  Xi{y)Xi{y)  +  C'm  Xm(2/)Xm(^)  +  C'i  {xf{y)Xm{y)  +  Xm(2/)Xf(^))  • 


(1.3) 


(1.4) 


We  note  that  (1.3)-(1.4)  give  the  pairwise  force  on  associated  with  a  unit  periodic  geom¬ 


etry.  In  summary,  the  function  a  in  (1.4)  is  given  by 


Cf,  if  y  and  y  are  in  the  fiber  phase 
o^iyii))  =  C'm,  if  y  and  y  are  in  the  matrix  phase 
Ci,  otherwise. 


In  equation  (1.3),  the  peridynamic  horizon  6  is  chosen  to  be  smaller  than  the  fiber  thickness 


in  the  unit  cell.  The  material  parameters  C{  and  Cm  are  intrinsic  to  each  phase  and  can  be 
determined  through  experiments.  Bonds  connecting  particles  in  the  different  materials  are 
characterized  by  C,  which  can  be  chosen  such  that  Cf  >  C  >  Cm  >  0,  see  p4]. 

The  microgeometry  associated  with  the  length  scale  e  is  obtained  by  rescaling  the  bond 
force  /short  as  follows.  For  x  E  fl, 

1  f X  x  +  ^\  ^  X 


/short(^,C,a:)  =  <  £2 


e 

0, 


otherwise. 


We  see  from  (1.4)  that  «(-,-)  is  given  by 


«  =  Ci  Xf{x)Xi{x)  +  Cm  Xm{x)xtn{x)  +  Ci  (Xf (a:)Xm(^)  +  xl,{x)x\{x)) ,  (1.5) 

where  xlix)  :=  Xf(f)  and  Xm{x)  :=  Xm(f)- 

The  peridynamic  equation  of  motion  for  this  model  is  given  by 


[X  —  X) 


dtu%x,t)=f 

Jh,s{x)^  \x-x\ 

supplemented  with  initial  conditions 


[X  —  X 


-(pY{x,  t)  —  u^{x,  t))  dx  +  h  4,  “5  4 


u^{x^  0)  =  (^x,  -  j  , 

dtu^{x,0)  =  • 


(1.7) 

(1.8) 


In  what  follows,  we  will  denote  by  s  a  real  number  such  that  |  <  s  <  cx).  In  (1.6)-(1.8), 


b{x,y,t)  is  in  C([0,T];  L^(f2  x  Y)^)  and  H-periodic  in  y  and  u^{x,y)  and  v^{x,  y)  are  in 
X  y)^  and  H-periodic  in  y. 
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Figure  4:  Long-range  bonds  (horizon  7)  and  short-range  bonds  (horizon  e6). 


In  the  second  model,  the  following  long-range  pairwise  force  is  added  to  the  short-range 
pairwise  force  of  the  hrst  model  (see  Figure]^ 

/l<.ng(r),0=  I 

I  0,  otherwise, 

where  7  is  a  prescribed  peridynamic  horizon.  Here  A  is  a  real-valued  function  dehned  by 


71 


C^,  otherwise. 


(1.9) 


where  denotes  the  angle  between  ^  and  a  line  parallel  to  the  hber  direction,  with 
0  <  <  |.  The  constants  Cf^  and  are  macroscopic  parameters  determined  through 

experiments,  see  [21IID]. 

Now  the  peridynamic  equation  of  motion  associated  with  the  total  pairwise  force  is  given 
by 


d^u^{x,t)  = 


/  N  (a:  -  xj  ®  (X  - 

/  X{x  —  x) - — - -  iu  ixR)  —  u  ixR))  dx 

I  rv*  ^y»  O 

'  H-y(x)  H  “^1 


+ 


1  f  X  x\  {x  —  x)  ®  {x  —  x) 

^  a  ' 


\  e  e 


\x  —  x\ 


{u^{x,  t)  —  u^{x,  t))  dx  (1-10) 


+  &  (x,  pt)  , 


Bacim  Alali  and  Robert  Lipton 


supplemented  with  initial  conditions 


0)  =  j  , 

dtu^{x,0)  =  . 


(1.11) 

(1.12) 


Remark  1.  The  hrst  model  follows  from  the  second  model  on  setting  A  =  0.  Thus  in 
Sections  |2]j^  we  will  often  present  our  results  and  analysis  for  the  second  model  only. 

In  the  third  model,  the  pairwise  force  is  given  by 


fiVy,^y,y) 


Isy  I 

0,  otherwise, 


where  y  eY  and  h  is  a  prescribed  peridynamic  horizon,  and  is  given  by 


aL{y,y  +  Q 


C'f 

C'm  141  ’ 


if  y  and  2/  +  4  in  the  hber  phase, 
and  4  is  parallel  to  n, 
otherwise. 


Here  4  is  the  Dirac  delta  distribution  concentrated  at  a  line  parallel  to  n.  The  function  ul 
can  be  written  in  terms  of  Xf  as  follows 


«l(2/,  1/  +  4)  =  C'f  141 4(4)  xi{y)xi{y  +  4)  +  Cm  141  (i  -  4(4)  xi{y)xi{y  +  4))  •  (i-i3) 


We  note  that  in  equation  (1.13),  Xiiv)  =  Xf(l/  +  4)  because  y  and  2/  +  4  both  he  on  a  line 
parallel  to  the  hber  direction  n. 

The  the  pairwise  force  dehned  on  D  is  given  by 


where  is  dehned  by 


al{x,x  +  ^)^-^ri,  1^1  <5 


0, 


otherwise. 


ai(x,x  +  4  =  Ci  1^14(4x14)  +  £ Cm  14  (1-4(4x14))-  (i-i4) 

The  peridynamic  equation  of  motion  for  this  model  is  given  by 

d?u^{x,t)  =  [  a^{x,x)— — }  ^ " - ^-iu^ixR)  —  u^ixR))  dx  (1-15) 

Jhs{x)  4  - 

supplemented  with  initial  data 

m'^(x,  0)  =  (x) ,  (1-16) 

C?iM^(x,  0)  =  (x) .  (1-17) 


Here  the  initial  data  vP  and  are  in  Lp(Q)^  with  1  <  p  <  cx)  and  the  loading  force  in 


equation  (1.15)  is  zero. 
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2  Multiscale  Analysis  and  the  Numerical  Scheme  for 
the  Short-Range  and  Long-Range  Bond  Model 


In  this  section,  we  present  the  multiscale  analysis  method  for  computing  the  deformation 
of  fiber-reinforced  composites  modeled  by  the  peridynamic  formulation.  This  is  done  for 
the  Short-Range  and  Long-Range  Bond  model  described  in  Section  1.2  The  method  de¬ 


livers  a  computationally  inexpensive  multiscale  numerical  scheme  for  the  analysis  of  these 
peridynamic  models  of  fiber-reinforced  materials.  It  consists  of  the  following  three  steps. 


1.  Macroscopic  Equation 

Compute  the  macroscopic  or  average  displacement  field  u^{x,t)  by  solving  a  peridy¬ 
namic  macroscopic  equation. 


2.  Cell— Problem 

Compute  the  micro-level  displacement  field  r{y,t)  by  solving  a  peridynamic  problem 
on  a  single  period  cell. 


3.  Downscaling 

The  displacement  field  of  the  oscillatory  peridynamic  equation  is  given  approximately 
by  superimposing  the  rescaled  micro-level  mechanical  responses  over  the  average  dis¬ 
placement  field,  i.e.,  u^Approx  —  +  r{x/e,t).  The  error  in  this  approximation  is 

shown  to  converge  in  norm  to  zero,  i.e.,  \\u'^{x,t)  —  UApprox{^y'^)\\  — >  0  as  £  — 0. 

In  the  following  subsections,  we  consider  four  cases  of  initial  and  loading  conditions. 
For  each  case,  we  present  the  macroscopic  equation,  the  cell-problem,  and  the  associated 
approximation.  The  results  provided  in  this  section  are  justified  in  Section 

For  convenience,  we  introduce  the  following  notation  for  the  average  of  a  periodic  func¬ 
tion.  Let  a  function  of  the  form  p{y),  p{x,y),  or  p{x,y,t)  be  R-periodic  in  the  variable  y. 
Its  average  over  Y  is  denoted  by 

p  = 

p{x)  = 
p{x,t)  = 

respectively.  For  future  reference,  we  let 

K=  X{x  —  x)— — ^  — —  dx  .  (2.1) 

JH-,ix)  \X  -  X\-^ 

By  the  change  of  variables  .^  =  a;  —  x,  it  is  easy  to  see  that  iF  is  a  constant  matrix,  which 
depends  on  the  macroscopic  parameters  7,  and 

For  future  reference,  we  will  adopt  the  notation  for  the  space  of  Lebesgue  p- 

integrable  functions  which  are  R-periodic.  Similarly,  Cper(R)  denotes  the  space  of  continuous 
R-periodic  functions.  Also  we  denote  by  C'°’^(r2)  the  space  of  Holder  continuous  functions 
with  exponent  /3,  where  0  <  /3  <  1. 


'Y 


'Y 


P{y)  dy, 
p{.x,y)dy,  or 
p{x,y,t)  dy, 


>Y 


(x  —  xd  6d  (x  —  X.) 
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2.1  First  Case 

In  this  section,  the  loading  force  and  initial  data  are  given  by 

b(^x,—,t^  =  l{x,t)  +  R  ,  (2.2a) 

=  uo(x)  +  Ml  j  ,  (2.2b) 

^  j  =  mo(2:)  +  ^^1  ,  (2.2c) 

where  I  G  C([0,T];  L®(t2)^),  R  is  in  with  .R  =  0,  Mq  and  Mq  are  in  U{VlY,  and  Wi 

and  Ml  are  in  L^^^iVY  with  -ui  =  hi  =  0.  Here,  12(f)  can  be  interpreted  as  a  residual  force. 
For  example,  such  forces  can  arise  from  the  differences  in  thermal  expansion  between  the 
two  materials. 


2.1.1  The  Macroscopic  Equation 

The  macroscopic  or  homogenized  peridynamic  equation  is  given  by 
dtu^{x,t)  = 


X{x  —  x) 


{x  —  x)®{x  —  x)jj  u 


\x  —  xP 


I  Hj{x) 

supplemented  with  initial  data 

M^(x,0)  =  uo{x),  dtu^{x,0)  =  vo{x). 


{u^{xR)  —  u^{x,t))  dx  +  l{x,t),  (2.3) 


(2.4) 


Here  the  macroscopic  displacement  is  the  weak  limit  of  the  sequence  of  displacements  u^. 
This  is  described  by  the  following  theorem. 


Theorem  2.1.  Let  he  the  solution  of  (1.10)-(1.12),  where  b,  ,  and  are  given  by 


(2.2).  Then  as  e  — >  0 


u^{xR)  — ^  u^{x,t)  weakly  in  L®(H  x  (0,T))^, 


where  E  C^([0,T];  L®(H)^)  is  the  unique  solution  of  (2.3)-(2.4). 

Moreover,  assume  that  I  E  C([0,T];  (7(12)^),  and  uq  and  vq  are  in  (7(12)^.  Then  is  in 
C2([0,T];C(H)3). 


2.1.2  The  Cell— Problem 

The  cell-problem  or  the  micro-level  peridynamic  equation  is  given  by 

=  [ 


I  Av-v)®{v-v)  I  I I 

- Y[y,t)-  r[y,t))  dy 


JHsiy) 

-K  r{y,t)  +  R{y), 
supplemented  with  initial  conditions 

r(y,0)  =  Ui(y),  d,r(y,0)  =  Vi(y). 


(2.5) 

(2.6) 


The  matrix  K  is  given  by  (2.1). 
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2.1.3  Downscaling 

The  macroscopic  displacement  together  with  the  rescaled  solution  of  the  cell  problem 
provide  the  approximation  to  the  actual  solution  given  by  u^Approx  —  +  r{x/e,t). 

This  is  expressed  in  the  following  theorem. 


Theorem  2.2.  Let  he  the  solution  of  (1.10)-(1.12),  where  b,  u^,  and  are  given  by 
(2.2).  Assume  that  I  E  C{[0,T];  0(0.)^),  and  Uq  and  Vq  are  in  C .  Then  for  almost  every 

t  e  (o,T), 


lim 

£— »'0 


u^{x,  t)  -  (^u^{x,  t)  +  r  t  j  j 


L'>(n)3 


=  0, 


(2.7) 


where  r  G  C'^([0,T];  LL^{Y)p  is  the  unique  solution  of  (2.5)-(2.6). 


Moreover,  assume  that  X  =  0  in  equation  (1.10).  Then,  for  i  G  (0,T)  and  uq,  vq,  and 
l{-,t)  in  the  error  in  (2.1)  is  estimated  by 


u^{x,t)  -  (u^{x,t)  +  r  <  Mi{t)e^, 

V  \e  //  inup 


L'>(0)3 

where  Mi{t)  is  independent  of  e.  The  function  Mi{t)  is  given  explicitly  in 
Section  \5.2.1[ 


(2.8) 


2.2  Second  Case 

In  this  section,  the  loading  force  and  initial  data  are  given  by 

6(a;,  pt)  =  F(^^,t^h{x), 


X 


where  F  G  C'([0,T];  L;,,(F)3x3)  and  h  G  L^(fi)= 


2.2.1  The  Macroscopic  Equation 

The  macroscopic  peridynamic  equation  is  given  by 
d'(u^{x,  t)  = 


(2.9a) 

(2.9b) 

(2.9c) 


I  H..,{x) 


X{x  —  x)— — — —  {u^{x,  t)  —  u^{x,  t))  dx  +  F{t)h{x),  (2.10) 


X  —  X 


supplemented  with  initial  data 

u^{x,0)  =  0,  dtu^{x,0)=0.  (2.11) 

Here  the  macroscopic  displacement  is  the  weak  limit  of  the  sequence  of  displacements  u^. 
This  is  described  by  the  following  theorem. 

Theorem  2.3.  Let  he  the  solution  of  (1.10)-(1.12),  where  b,  ,  and  are  given  by 


(2.9).  Then  as  £  — >  0 


u^(x,t)  — >  u^{x,t)  weakly  in  L^{Vt  x  (0,T))  , 


where  G  C^([0,T];  L®(12)^)  is  the  unique  solution  of  (2.10)-(2.11). 
Moreover,  assume  that  h  G  C'(r2)^.  Then  is  in  C^([0,T];  (7(12)^). 


Bacim  Alali  and  Robert  Lipton 


12 


2.2.2  The  Cell— Problem 

The  micro-level  peridynamics  is  given  by  the  following  equations.  For  j  =  1,  2,  3, 
dy{y,t)  =  [  a{y,y)^^  ^  (r^iyR)  -  r^ivR))  dy 

jHsiy)  \[y-y)\ 

-Kr^{y,t)  +  {F^{y,t)-F^{t)),  (2.12) 


supplemented  with  initial  conditions 

r^(l/,0)  =  0,  dtF{y,0)=0. 


(2.13) 


In  (2.12),  F^{y,t)  and  F^{t)  denote  the  columns  of  the  matrices  F{y,t)  and  F{t),  respec¬ 
tively.  The  matrix  K  is  given  by  (2.1). 


2.2.3  Downscaling 

The  macroscopic  displacement  together  with  the  rescaled  solution  of  the  cell  problem 
provide  an  approximation  to  the  actual  solution  .  This  is  expressed  in  the  following 
theorem. 


Theorem  2.4.  Let  he  the  solution  of  where  h,  ,  and  are  given  by 

(2.9).  Assume  that  h  G  (7(0)^.  Then  for  almost  every  t  G  (0,T), 


lim 


u^{x,  t)  —  (  u^(x,  t)  -|-  F  ,  t )  hj{x) 


X 


=  0, 


Lqr2)3 


where  F  G  ^^([OjT];  is  the  unique  solution  of  (2.12)-(2.F^. 

Moreover,  assume  that  X  =  0  in  equation  (1.10).  Then,  for  t  G  (0,  T)  and 
h  G  the  error  in  (2. If)  is  estimated  by 


F{x,  t)  —  (  u^{x,  t)  +  F  t )  hj{x) 


X 


<  M2(t)£^ 


L^inp 


where  M2(t)  is  independent  of  e.  The  function  M2(t)  is  given  explicitly  in 
Section  \5.2.1\ 


(2.14) 


(2.15) 


2.3  Third  Case 

In  this  section,  the  loading  force  and  initial  data  are  given  by 


=  0, 

(2.16a) 

“”(-f) 

=  F(^^^h{x), 

(2.16b) 

(..  1) 

=  0, 

(2.16c) 

where  F  G  Lpg^(y)^^^  and  h  G  L^iVt)'^. 
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2.3.1  The  Macroscopic  Equation 

The  macroscopic  peridynamic  equation  is  given  by 

dfu^{x,t)=  [ 


supplemented  with  initial  data 

{x,0)  =  Fh{x),  dtu^{x,0)  =  0. 


(2.17) 


(2.18) 


Here  the  macroscopic  displacement  is  the  weak  limit  of  the  sequence  of  displacements  u^. 
This  is  described  by  the  following  theorem. 


Theorem  2.5.  Let  he  the  solution  of  (1.10)-(1.12),  where  b,  ,  and  are  given  by 
(2.16).  Then  as  £  — >•  0 


u^{xR)  — 1-  u^{x,t)  weakly  in  L^{Ll  x  (0,T))^, 


where  G  L^{Vt)‘^)  is  the  unique  solution  of  (2.11)-(2.18). 

Moreover,  assume  that  h  G  Then  is  in  C^([0,T];  C(r2)^) 


2.3.2  The  Cell— Problem 

The  micro-level  peridynamics  is  given  by  the  following  equations.  For  j  =  1,  2,  3, 


dy{y,t)  = 


a{y,y)^^  dy 


JHsiy) 

-K  F{y,t), 

supplemented  with  initial  conditions 

F{y,0)  =  F^{y)-F^,  dtF{y,0)  =  0 


(2.19) 


(2.20) 


In  (2.20),  F^{y)  and  F^  denote  the  columns  of  the  matrices  F{y)  and  F,  respectively. 
The  matrix  K  is  given  by  (2.1). 


2.3.3  Downscaling 

The  macroscopic  displacement  together  with  the  rescaled  solution  of  the  cell  problem 
provide  an  approximation  to  the  actual  solution  .  This  is  expressed  in  the  following 
theorem. 


Theorem  2.6.  Let  be  the  solution  of  (1.10)-(1.12),  where  b,  ,  and  are  given  by 
(2.16).  Assume  that  h  G  (7(0)^.  Then  for  almost  every  t  G  (0,T), 


lim 


u^{x,  t) 


(2.21) 


Bacim  Alali  and  Robert  Lipton 


14 


where  E  C^([0,T];  Lpg^(y)^)  is  the  unique  solution  of  (2.19)-(2.2(^. 

Moreover,  assume  that  X  =  0  in  equation  (1.10).  Then,  for  t  E  (0,  T)  and  h  E 


the  error  in  (2.21)  is  estimated  by 

u^{x,  t)  —  (u^{x,  t)  +  r^  hj{x) 


i=i 


<  Ms{t)e^, 


L‘>(0)3 


where  M^{t)  is  independent  of  e.  The  function  M^it)  is  given  explicitly  in 
Section  \5.2.tA 

2.4  Fourth  Case 

In  this  section,  the  loading  force  and  initial  data  are  given  by 


where  F  E  and  h  E 

2.4.1  The  Macroscopic  Equation 

The  macroscopic  peridynamic  equation  is  given  by 


dfu  {x,t)  =  /  X{x  —  x) 
Jhj{x) 


(x  -  x)  (gi  (x  -  x)  jj 


\x  —  xP 


{u  {x,  t)  —  u  {x,  t))  dx, 


supplemented  with  initial  data 

u^{x,0)  =  0,  dtu^  {x,0)  =  Fh{x). 


(2.22) 


=  0, 

(2.23a) 

“”(-f) 

=  0, 

(2.23b) 

=  F  (1)  h{x), 

(2.23c) 

(2.24) 


(2.25) 


Here  the  macroscopic  displacement  is  the  weak  limit  of  the  sequence  of  displacements  u^. 
This  is  described  by  the  following  theorem. 


Theorem  2.7.  Let  be  the  solution  of  (1.10)-(1.12),  where  b,  m,  and  are  given  by 


(2.23).  Then  as  e  ^  0 


u^{x,t)  — >•  u^{x,t)  weakly  in  L^{Ll  x  (0,T))'*, 


where  E  C^([0,T];  is  the  unique  solution  of  (2.24)-(2.25). 

Moreover,  assume  that  h  E  C{Cl)^.  Then  is  in  C^([0,T];  (7(12)^). 
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2.4.2  The  Cell— Problem 

The  micro-level  peridynamics  is  given  by  the  following  equations.  For  j  =  1,  2,  3, 


dy{y,t)  =  [  a{y,y)^^  ^  {r^{y,t)  -  r\y,t))  dy 

jHsiy)  \[y-y)\ 

-K  r^{y,t), 

supplemented  with  initial  conditions 

{y,  0)  =  0,  dtr^  (l/,  0)  =  (y)  -  P . 


(2.26) 


(2.27) 


In  (2.27),  F^{y)  and  denote  the  columns  of  the  matrices  F{y)  and  F,  respectively. 
The  matrix  K  is  given  by  (2.1). 


2.4.3  Downscaling 

The  macroscopic  displacement  together  with  the  rescaled  solution  of  the  cell  problem 
provide  an  approximation  to  the  actual  solution  .  This  is  expressed  in  the  following 
theorem. 


Theorem  2.8.  Let  he  the  solution  of  (1.10)-(1.12),  where  h,  P,  and  are  given  by 
(2.23).  Assume  that  h  G  (7(11)^.  Then  for  almost  every  t  G  (0,T), 


lim 


P{x,  t)  —  (  u^{x,  t)  +  F  ,  t )  hj{x) 

i=i 


X 


=  0, 


(2.28) 


Lyo)3 


where  F  G  C^([0,T];  L^pY)^)  is  the  unique  solution  of  (2.26)-(2.21^. 

Moreover,  assume  that  X  =  0  in  equation  (1.10).  Then,  for  t  G  (0,  T)  and  h  G 
the  error  in  (2.28)  is  estimated  by 


P{x,  t)  —  (  u^{x,  t)  +  F  t )  hj{x) 


X 


<  Mp)e^, 


(2.29) 


where  Mp)  is  independent  of  e.  The  function  Mp)  is  given  explicitly  in 
Section  15.  ^.51 


3  Existence  and  Uniqueness  Results  for  the  Peridy- 
namic  Equation 


In  this  section,  we  make  use  of  semigroup  theory  of  operators  to  study  the  existence  and 
uniqueness  of  ( 1.10 )-( 1.12 ).  We  begin  by  introducing  the  following  operators.  For  v  G  L^pL)^, 
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with  I  <  s  <  cx),  let 

Al,iv{x) 

Al,2v{x) 

Also  we  set 


Ax-x  ®[x-x 
A{x  —  X) - - v{x)  ax, 


\x  —  x\ 


Ihax) 


Ax-x  ®[x-x 
A{x  —  X) - - ax  v[x), 


\x  —  x\ 


1  f  X  x\  {x  —  x)  {x  —  x) 


'  H^s{x) 


\X  —  X\ 


1  f  X  x\  {x  —  x)  ^  {x  —  x) 


'  H^s{x) 


\x  —  x\ 


v(x)  dx, 


dx  v{x). 


Al 

A^s 

A" 


Al,1  —  Al^2 

^S,l  ^S,2 

Al  +  Ag. 


5 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


(3.5) 

(3.6) 

(3.7) 


Then  by  making  the  identihcations  u^{t)  =  u‘^{-,t)  and  =  b{-,  -,t),  we  can  write  (1.10)- 
(1.12)  as  an  operator  equation  in  L^{Q)^ 


(  vA{t)  =  A^u^{t)  AhF{t),  tG[0,T] 

I  M^(0)  =  ul, 

y  -U^(O)  =  Ug. 

or  equivalently,  as  an  inhomogeneous  Abstract  Cauchy  Problem  in 

uiSif  X  uisif 

I  =  A^U^{t)  +  te[0,T] 

I  Cl'(O)  =  f/„'. 

where 

=  (  “:j(J  )  .  l/o'  =  (  :?  )  .  B^t)  =  (  )  ,  and  A« 


(3.8) 


(3.9) 


0 

A” 


Here  I  denotes  the  identity  map  in  L®(r2)^. 

Proposition  3.1.  Let  |  <  s  <  cxd  and  assume  that  ¥  G  ^([OjT];  L^{VtY).  Then 

(a)  The  operators  and  are  linear  and  hounded  on  L^{Ll)^  and 

L^{Q)^  X  respectively.  Moreover,  the  bounds  are  uniform  in  e. 

(b)  Equation  (3.f^  has  a  unique  classical  solution  in 

C^([0,T];  L®(r2)^  X  which  is  given  by 


U%t)=e^^'U^+  /  te[0,T], 


(3.10) 
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where 


e  =  >  — 
n\ 

n=0 


_  ~  4-fl 

E  (A')" 


(3.11) 


Moreover,  equation  (3.8)  has  a  unique  classical  solution  G  C^([0,T];  which 

is  given  by 


u^{t)  =  cosh  + '/(M  sinh 

j  Anh.  {t)  dr 


£  UG 


with  the  notation 


cosh  {t'/^J  := 


°°  ^2n 


n=0 


(2n) 


sinh  {t'/(M'\  := 


^2nH-l 

^  (2n  +  1)! 

n=0  ^  ^ 


(A" 


(3.12a) 

(3.12b) 

(3.12c) 


The  sequences  {u^)e>o,  {id)e>ci,  and  {u^)e>o  are  hounded 
in  h“([0,T];  L^(fi)3). 

Proof.  Part  (a).  It  is  clear  that  the  operators  A|;^,  A|2;  and  Al^2  are  linear.  So 

we  begin  the  proof  by  showing  that  A|^  and  A52  are  uniformly  bounded  on  L^(f2)^  for 
I  <  s  <  cxo.  Let  V  G  L^(fl)^.  Then  by  the  change  of  variables  x  =  x  +  ez  in  (3.3)  we  obtain 


^l,A(a;)  = 


'Hs{0) 


X  X 


a  +  z 


e  e 


z  ®  z 


v{x  +  ez)  dz. 


(3.13) 


Let  ttmax  =  max  a{y,y').  Then  by  taking  the  Euclidean  norm  in  (3.13),  we  see  that 

y,y'eY 


|^|,A(a;)|  <  CTmax 

E  C^max 


Hs(0) 


— j-  \v{x  +  ez)  \  dz 


'Hs{0) 


-7  dz 

c' 


1/s’ 


'Hs{0) 


\v{x  +  ez)\^  dz 


l/s 


(3.14) 


where  Holder’s  inequality  was  used  in  the  second  inequality,  with  l/s  +  l/s'  =  l 

and  1  <  s'  <  3.  By  changing  the  variable  of  integration  back  to  x  in  the  second  integral, 

and  then  taking  the  limit  as  e  — >■  0,  we  see  that 


'Hs{0) 


\v{x  +  ez)f  dz  =  “§■  / 

^  J  Hssix) 


|n(a;)|^  dx 

— >•  \Hs{x)\  |f(a;)|^,  a.e.  x,  (3.15) 

where  we  have  used  Lebesgue’s  Differentiation  Theorem  to  evaluate  this  limit.  On  the  other 
hand,  we  observe  that  the  hrst  integral  in  (3.14)  is  hnite  because  s'  <  3.  Therefore,  it  follows 
from  (3.14)  and  (3.15)  that 

|41|^A(x)|  <  Mi|r;(a;)|, 
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for  some  real  number  Mi  >  0  which  is  independent  of  e.  It  follows  that 

<  Ml  ||n||Ls(n)3, 

which  shows  that  the  operator  Ag  ^  is  uniformly  bounded.  Similarly,  A^g  2  can  be  written  as 

z  z 


= 


/  fx  X 
/  a  +  z 

'Hs(0)  £ 


dz  v{x). 


(3.16) 


Thus 


'Hs{0) 


—  dz  |n()r)|, 


from  which  the  boundedness  of  A^g2  immediately  follows.  Combining  these  results  shows 
that  y4|,  which  is  given  by  A^g^  —  A521  ^  uniformly  bounded  operator  on 

Next  we  show  that  the  linear  operator  Ai  =  A^^i  —  Al^2  is  bounded  on  Let 

Amav  =  max  A(^).  Then  by  taking  the  Euclidean  norm  in  (3.1),  we  see  that 


\Al,iV{x)\  <  Ar 


<  A, 


'  H^{x) 


\X  —  X\ 


|n(a;)|  dx 


I/s' 


l/s 


\x  —  x\ 


dx 


f  H-y{x) 


|n(a;)|*  dx 


(3.17) 


where  Holder’s  inequality  was  used  in  the  second  inequality,  with  l/s  +  l/s'  =  l 

and  1  <  s'  <  3.  By  the  change  of  variables  ^  =  L  —  x,  it  is  easy  to  see  that  the  hrst  integral 

in  (3.17)  is  independent  of  x  and  hnite  because  s'  <  3.  Therefore  from  (3.17)  we  obtain 

I/s' 


II^L, 1^11^5(0)3  <  Ar 


dz 


TllLqo)3- 


This  shows  that  Ai^i  is  bounded  on  The  boundedness  of  Al^2,  which  is  given  by 

(3.2),  is  clear.  Therefore  A^  is  bounded  on 


Since  A^  =  Al  +  Ag,  we  conclude  that 


(3.18) 


for  some  real  number  M  >  0  which  is  independent  of  e. 

The  operator  is  clearly  linear,  thus  it  remains  to  show  that  this  operator  is  uniformly 
bounded  on  x  To  see  this,  we  let  {v,w)  G  x  The  norm  in  this 

Banach  space  is  given  by 


(n,  ■?n)||L5(o)3xL5(o)3  —  11^11^5(0)3  +  ||rc||L5(o)3. 


We  note  that 


A" 


w 


0  / 
0 


V 

w 


w 

A^v 
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Thus  by  taking  the  norm,  we  obtain 


||A'^(u,  r(;)||Ls(n)3xL'>(o)3  —  ||'M^||L«(n)3  + 

<  ||'U^||L«(r2)3  +  ||A^||  ||'y||L^(n)3- 


From  (3.19)  and  since  we  may  assume  that  M  >  1  in  (3.18),  it  follows  that 


||A^(u, r(;)||is(Q)3xL'’(r2)3  <  Af||(u,r(;)||is(n)3xL'*(o)3, 


(3.19) 


(3.20) 


completing  the  argument. 

Part  (b).  We  have  seen  from  Part  (a)  that  is  a  bounded  linear  operator  on  the  Banach 
space  L^(r2)^  x  L^(r2)^.  Also,  since  ¥  is  in  C([0,T];  by  assumption,  it  follows  that 

=  (0,6^)  is  in  C([0,T];  x  From  these  facts,  it  follows  from  the  theory  of 

semigroups  thatQ 

1.  The  operator  A*^  generates  a  uniformly  continuous  semigroup  on  L^(n)3  x 

where  is  given  by  (3.11). 


2.  The  inhomogeneous  Abstract  Cauchy  Problem  (3.9)  has  a  unique  classical  solution 


e  Cn[0,T];  X  L%nf)  which  is  given  by  (3.10). 


It  immediately  follows  from  (|^  that  the  second  order  inhomogeneous  Abstract  Cauchy  Prob¬ 
lem  (3.8)  has  a  unique  classical  solution  e  C^([0,T];  It  remains  to  show  that 


is  given  explicitly  by  (3.12).  To  see  this,  we  begin  by  the  following  observations  which  can 


be  easily  shown  using  mathematical  induction.  For  n  =  0, 1, 2, . . .,  we  have 


0  I 
0 


2n 


(A 


e\n  Q 

{A^y 


0 


(3.21) 


0  I 
A^  0 


2n+l 


0  (A^)^ 

(^£)n+l  Q 


From  (3.11)  and  by  using  these  two  equations  we  see  that 

( 


j-hr 


n=0 


¥  f  0  I 

n!  V  0 


^  ^2n 


E 


(2n)! 

n=0  ^  ’ 

°°  ^2n+l 


(W)^ 


(W) 


“  j-2n+l  \ 

n=0 

^  f2n 

(^£)n 


E 

n=0 


(2n)! 


Equation  (3.12)  follows  from  equations  (3.10)  and  (3.23),  and  the  fact  that 


= 


u 

if 


(3.22) 


(3.23) 


'^see  for  example  unma. 
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Part  fc).  We  recall  that 


Uq{x)  :=  -) 

X, 

or'  - 


Vnix)  :=  v^{x,  — 


Also  by  assumption  u^{x,y),  v^{x,y)  are  in  Lp^^(Y)^).  Therefore  we  see  that 


|'Wo||i,=  (o)3  <  ||M°|U=(0;Ls^^(y)3)  ■—[  I  y  1“ 


kolUhr®  -  lk°lk«rTierr)3)  •- 


n  JY 


\  l/« 

^{x,y)\''  dydxj  , 
\  l/s 

v^{x,y)Y  dydx\ 


Thus  Uq  and  Ug  are  uniformly  bounded  in  L^iVL)^,  which  implies  that  Uq  is  uniformly 
bounded  in  L^{Q)^  x  Similarly  we  can  show  that  for  t  G  [0,T],  is  uniformly 

bounded  in  L®(r2)^.  Since  R{t)  is  continuous  in  t,  it  follows  that  ¥  is  uniformly  bounded  in 
r([0,T];  which  implies  that  is  uniformly  bounded  inC'([0,T];  (Q,)^  x  . 

Next  we  note  that 


<  B 


<  (3.24) 

where  in  the  last  inequality  we  have  used  the  fact  that  M  is  uniformly  bounded.  Taking  the 


norm  in  both  sides  of  (3.10)  and  by  using  (3.24),  we  obtain 

rt 


(3.25) 


for  some  positive  numbers  Mi,  M2,  and  M.  This  implies  that  is  uniformly  bounded 
in  L°°([0,T];  x  L^(r2)^).  Therefore  the  sequences  {u^)e>o  and  ('U^)£>o  are  bounded 

in  L“([0,T];  Finally,  it  follows  from  equation  (3.8)  that  the  sequence  ('ii'^)e>o  is 

bounded  in  L°°([0,T];  L^(r2)^),  completing  the  proof.  □ 


4  Two-Scale  Convergence  and  the  Two-Scale  Limit  Equa¬ 
tion 


The  aim  of  this  section  is  to  identify  the  two-scale  limit  of  the  peridynamic  initial-value 
problem  ( 1.10 )-( 1.12). 


4.1  Two-Scale  Convergence 

We  begin  by  dehning  two-scale  convergence  and  recalling  some  results  from  two-scale  con¬ 
vergence.  In  the  subsequent  discussion,  we  will  often  refer  to  the  following  function  spaces 

/C  =  {'0  e  ^“(M^  X  F),  il){x^y)  is  F-periodic  in  y}, 

J  =  {ip  &  ^“(M^  X  F  X  M+),  ip{x,y,t)  is  F-periodic  in  y}, 

Q  =  [w  E  C'^([0,T];  L^{Q  X  F)^),  w{x,y,t)  is  F-periodic  in  y,  and  3/2  <  s  <  cx)}. 
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Let  p  and  p'  be  two  real  numbers  such  that  1  <  p  <  cx)  and  1/p  +  1/p'  =  1. 

Definition  4.1  (Two-scale  convergence  ||18L  IT]).  A  sequence  (n^)  of  functions  in  L^{Q), 
is  said  to  two-scale  converge  to  a  limit  v  G  WiVL  x  Y)  if,  as  £  — 0 


v‘^ (x , —  j  dx  — >  /  v{x,y)'i/{x,y)  dxdy 

'  -'nxY 


(4.1) 


'  2 

for  all  f)  E  (12;  Cper{Y)).  We  will  often  use  ^  v  to  denote  that  (v^)  two-scale  converges 
to  V. 


If  the  sequence  (v^)  is  bounded  in  L'P{Vt)  then  L"^' [Vt]  Cper{Y))  can  be  replaced  by  /C  in 
Definition  (4.1)  (see  [T^L 

The  following  are  well-known  results  on  two-scale  convergence,  which  can  be  found  in 

[19]. 


Proposition  4.2.  If  {v^)  converges  to  v  in  LP(f2)  then  (n^)  two-scale  converges  to  v{x,y)  = 
v{x). 

Proposition  4.3.  If  fj  E  K,  then  'ip{x,  |)  two-scale  converges  to  'i/>{x,y). 

Proposition  4.4.  Let  {v^)  be  a  sequence  in  LP(f2)  which  two-scale  converges  to  v  E  L^(12  x 
Y).  Then 


v^{x)'ip  (x,  -j  dx  — >•  /  v{x,y)i/>{x,y)  dxdy, 

’  -'nxY 


for  every  x/  of  the  form  xjj{x,y)  =  xjji{x)xjj2{y) ,  where  xpi  E  (H)  and  xjj2  G  L^f/fY),  with 
1  <  r  <  oo  and  l/r-|-l/r'  =  l. 

Proposition  4.5.  Let  {v^)  be  a  sequence  in  L^(f2)  which  two-scale  converges  to  v  E  L^(12  x 
Y).  Then  as  e  — >•  0 

y  v{x,y)dy  weakly  in  L^{Il). 


Definition  4.1  is  motivated  by  the  following  compactness  result  of  Nguetseng,  see 

Theorem  4.6.  Let  (v^)  be  a  bounded  sequence  in  L^{fl).  Then  there  exists  a  subsequence 
and  a  function  v  E  L^(12  x  Y)  such  that  the  subsequence  two-scale  converges  to  v. 


For  the  time-dependent  problems  studied  in  this  work,  we  slightly  modify  the  above  two- 
scale  convergence  definition  and  results  to  allow  for  homogenization  with  a  parameter,  see 
EE]-  Here  the  parameter  is  denoted  by  t. 

Definition  4.7.  A  sequence  {v^)  of  functions  in  L^{VL  x  (0,T)),  is  said  to  two-scale  converge 
to  a  limit  v  E  L^(12  x  T  x  (0,  T))  if,  as  e  —>■  0 


'Qx(0,T) 


v‘^{x,t)xp  ^x,—,t^  dxdt 


>nxYx{0,T) 


v{x,y,t)xp{x,y,t)  dxdydt  (4.2) 


for  all  xp  E  J . 
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Theorem  4.8.  Let  be  a  bounded  sequence  in  x  (0,T)).  Then  there  exists  a  subse¬ 
quence  and  a  function  v  G  xYx  (0,T))  such  that  the  subsequence  two-scale  converges 
to  V. 


The  proof  of  this  result  is  essentially  the  same  as  the  proof  of  Theorem  4.6  A  slight 


variation  of  Theorem  4.8  can  be  found  in  |Hj  and  [^. 

The  following  is  a  direct  consequence  of  Dehnition  |4.7  and  the  dehnition  of  weak  conver¬ 
gence. 


Proposition  4.9.  Let  (n^)  be  a  bounded  sequence  in  L^{VL  x  (0,T))  that  two-scale  converges 
to  V  &  LP{Ll  xY  X  (0,T)).  Then  as  e  — >•  0 


>  y  v{x,y,t)dy  weakly  in  L'^{Q  x  {0,T)). 

Finally,  we  state  the  following  well-known  result  on  the  weak  limit  of  oscillatory  periodic 
functions,  which  can  be  found  in  [B]. 

Proposition  4.10.  Let  h  G  L^{VL)  be  a  Y -periodic  function,  where  1  <  q  <  oo.  Set  h‘^{x)  = 
h(|)  for  X  e  Ll.  Then  as  e  ^  0, 


if  ^  Y  q  <  oo,  and 
if  q  =  oo. 


h^ 


weakly  in 


weakly-*  in 


(4.3) 

(4.4) 


4.2  The  Two-Scale  Limit  Equation 


In  this  section,  we  use  two-scale  convergence  to  identify  the  limit  of  ( 

1.10 

)-( 

1.12 

.  We  observe 

that  the  loading  force  and  initial  data  given  by  equations  ( 

2.2 

),  ( 

2.9) 

(2.16 

),  or 

2.23 

),  satisfy 

the  following 


j  v^{x,y). 


(4.5a) 

(4.5b) 

(4.5c) 


We  note  that  from  Proposition  3.1[c)  and  Theorem  4.8  it  follows  that,  up  to  some  subse¬ 
quences,  u,  if  ^  u*,  and  il^  u**,  where  u,  u*,  and  u**  are  in  L'^([0,T];  L^(f2  x  Y)^). 
We  shall  see  later  that  u{x,y,t)  is  uniquely  determined  by  an  initial  value  problem.  There¬ 
fore  u  is  independent  of  the  subsequence,  and  the  whole  sequence  (n^)  two-scale  converges 
to  u. 
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In  order  to  identify  the  two-scale  limit  of  (1.10),  we  mnltiply  both  sides  by  a  test  fnnction 


■^(x,  -,t),  where  ip{x,y,t)  is  F-periodic  in  y  and  is  snch  that 
^  g  X  y  X  M)^,  and  integrate  on  x  M+ 


'  OxM+ 


'  OxE+ 


d^u^(x,  t)  ■  Ip  j  dxdt 

(Al  +  Ag)u^{x,  t)  +  b  j  ■  V'  (^2;,  dxdt 

After  integrating  by  parts  twice,  we  obtain 

/  u^{x,t)  ■  d'f'ip  (x, —,tj  dxdt—  /  dtu^{x,0)  ■  pj  (x,  —,0j  dx 

JnxR+  \  e  /  Jq  \  e  / 

+  j  u^{x,t))  ■  dtp)  ^x,  — ,oj  dx 

{Al  +  Ag)u^{x,  t)  +  b  ^x,  dxdt 


J  Q,xR+ 

By  letting  e  — >■  0  we  obtain 


/  u{x,y,t)  ■  d‘lp){x,y,t)  dxdydt  —  /  v^{x,y)  ■  p){x,y,b)  dxdy 

IO.xYxR+  JqxY 


InxY 


=  lim 

£— ^0 


M°(x,  y)  ■  dtp){x,  y,  0)  dxdy 
/  {Al  +  Ag)u^{x,  t)  ■  p)  (x,  — ,  tj  dxdt 
b{x,  y,  t)  ■  p){x,  y,  t)  dxdydt 


lnxYx'R+ 


(4.6) 


For  i  =  1,  2,  3,  we  extend  Ui{x,y,t)  by  periodicity  from  fl  x  F  x  (0,T)  to 

fl  X  X  (0,T).  We  will  use  the  following  lemma  to  compute  the  limit  on  the  right  hand 


side  of  (4.6) 


Lemma  4.11.  Let  w  be  in  Lpg^(F)^)  and  define 

(x  —  x)  ®  (x  —  x) 


BLw{x,y) 

=  f  \{x  — . 

J  Hj{x) 

Bsw{x,y) 

=  /  a{y,y) 

JHsiy) 

Then  as  e  ^  0, 

(a)  ALu‘'{x,t) 

2 

^  BLu{x,y,t). 

\X  —  \J  Y 

{y-y)(»{y-  y) 


w(x,  y')  dy'  —  w{x,  y)  )  dx, 


\y-y\' 


{w{x,y)  -  w{x,y))  dy. 


Moreover,  the  operator  Bl  is  linear  and  hounded  on  L‘^{Vt]  Lp^fiY)^). 

(b)  A%ufix,t)  ^  Bsu{x,y,t). 

Moreover,  the  operator  Bs  is  linear  and  bounded  on 
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The  proof  of  this  lemma  is  provided  at  the  end  of  this  section. 

Using  Lemma  (4.11)  and  Lebesgue’s  dominated  convergence  theorem,  it  follows  that 

{Al  +  Ag)u^{x,  t)  ■  t  j  dxdt 


lim 

£— >0 


f!xM+ 


'Ox YxR+ 


{Bl  +  Bs)u{x,  y,  t)  ■  '(/’(x,  y,  t)  dxdydt. 


Thus  (4.6)  becomes 


'nxYxM.+ 


u{x,  y,  t)  ■  d‘lil){x^  y,  t)  dxdydt 


v^{x,y)  ■  'ijj{x,y,0)  dxdy 


>nxY 


M°(x,  y)  ■  dt'iplx,  y,  0)  dxdy 


inxY 


'OxYxE+ 


((Bi  +  Bs)u(x,y,t)  +  b{x,y,t))  ■■il}(x,y,t)dxdydt 


(4.7) 


We  shall  see  from  Lemma  4.13  provided  before  the  end  of  this  section,  that  u  has  two 
classical  partial  derivatives  with  respect  to  t,  for  almost  every  t,  and  the  initial  conditions 
supplementing  (4.7)  are  given  by 


'x(x,y,0)  =  dtu(x,y,0)  =  v°(x,y). 


(4.8) 


Thus  by  integrating  by  parts  twice,  equation  (4.7)  becomes 

d^u{x,  y,  t)  ■  'il^i^x,  y,  t)  dxdydt 


'QxYxR+ 


{{Bl  +  Bs)u{x,  y,  t)  +  6(x,  y,  t))  ■  '0(x,  y,  t)  dxdydt  (4.9) 

xY  X  M)^  for  which  'ijj{x,y,t)  is  U-periodic  in 


./Ox  Yx]R+ 

Since  this  is  true  for  any  function  'ip  ^  C\ 
y,  we  obtain  that  for  almost  every  x,  y,  and  t 

dfu(x,  y,  t)  =  Bu{x,  y,  t)  +  b(x,  y,  t), 


(4,10) 


where  B  =  Bl  +  Bs-  It  follows  from  Lemma  [4.11|  that  i?  is  a  bounded  linear  operator  on 
L®  (y)^).  Therefore,  the  initial  value  problem  given  by  (4.10)  and  (4.8),  interpreted 


as  a  second-order  inhomogeneous  abstract  Cauchy  problem  dehned  on  L®g^(y)^),  has 

a  unique  solution  u  E  Q. 

The  following  summarizes  the  results  of  this  section. 


Theorem  4.12.  Let  [yd]  he  the  sequence  of  solutions  of  Then 

2 

u  where  u  E  Q  is  the  unique  solution  of 

(x  —  x)  0  (x  —  x) 


d^u{x,y,t)  =  /  A(x-x)- 


'H-y{x) 


X  —  X 


m(x,  y' ,  t)  dy'  —  m(x,  y,  t)  )  dx 


'Y 


[  a{y,y)— — — —{u{x,yR)-u{x,yA))dy 

’Hsiy)  \y  -  yr 


(4.11) 


+  b{x,y,t), 


Bacim  Alali  and  Robert  Lipton 


25 


supplemented  with  initial  conditions 

u{x,y,0)  =  u°{x,y), 
dtu{x,y,0)  =  v^{x,y). 

Lemma  4.13.  Let  t  G  [0,T]  and  define 


/  u"(x,y,l)dldT +  tu’(x,y,0) +u{x,y,l)). 


(4.12) 

(4.13) 


(4.14) 


0  Jo 


Then  g  is  in  x  Y  x  (0,T))^,  twice  differentiable  with  respect  to  t  almost  everywhere, 
and  satisfies 

(a)  For  almost  every  x,y,  andt,  g{x,y,t)  =  u{x,y,t),  dtg{x,y,t)  =  u*{x,y,t), 
and  dlg{x,  y,  t)  =  u**  {x,  y,  t) . 

(b )  For  almost  every  x  and  y 

9{x,  y,  0)  =  u{x,  y,  0)  =  u^{x,  y), 
dtg{,x,y,b)  =  u*{x,y,b)  =  v^{x,y). 

Proof.  Part  (a).  Let  'i/ji{x,  y)  be  in  x  Y)^  and  F-periodic  in  y,  and  let  0  be  in  C'“(M+). 

Then  by  using  integration  by  parts,  we  see  that 


/f!xM+ 


dtU^{x,  t)  ■  tpi  (^x,  —  j  fit)  dxdt 


'  OxM+ 


t)  ■  01  j  (fit)  dxdt. 


2 

Sending  £  to  0  and  using  the  fact  that,  up  to  a  subsequence,  dtu^  u*,  we  obtain 


lnxYxR+ 


JqxYxR+ 

Since  this  holds  for  every  0i  we  conclude  that 


u*ix,  y,  t)  ■  01  {x,  y)  0(f)  dxdydt 
uix,  y,  t)  ■  01  (x,  y)  0(f)  dxdydt. 


u*{x,y,t)(j){t)  dt  =  —  /  u{x,y,t)f{t)  dt, 


(4.15) 


for  almost  every  x  and  y  and  for  every  0  G  C'“(M''').  Similarly,  by  using  the  fact  that,  up  to 
a  subsequence,  d'^u^  ^  u** ,  we  see  that 


Ir+ 


u**  [x,  y,  f )0(f )  dt=  uix,  y,  f )0(f )  dt, 


(4.16) 


for  almost  every  x  and  y  and  for  every  0  G  C'“(M+).  We  note  that  from  (4.14)  it  is  easy  to 
see  that  g  is  twice  differentiable  in  f  almost  everywhere  and  satishes 


dtgix,y,t)  =  /  u**ix,y,T)dT  +  u*ix,y,0), 

Jo 


(4.17) 

(4.18) 
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We  will  use  these  facts  together  with  (4.15)  and  (4.16)  to  show  that  dtg  =  u*  almost  every¬ 
where  and  g  =  u  almost  everywhere. 

For  0  e  we  have 


dtg{x,y,t)^{t)dt  = 


dtg{x,y,t)(p{t)dt 


=  -  u**{x,y,t)(f){t)dt 


=  —  u{x,y,t)(j){t)  dt 

Jr+ 


u*{x,  y,  t)0(t)  dt 


where  (4.18)  and  (4.16)  were  used  in  the  second  and  third  steps,  respectively.  Thus  we  obtain 

f  {dtg{x,y,t) -u*{x,y,t))^{t)dt  =  0,  (4.19) 


for  every  <p  G  C“(]R’'').  Since  dtg{x,y,0)  =  u*{x,y,0),  we  conclude  from  (4.19)  that 
dtg{x,y,t)  =  u*{x,y,t)  almost  everywhere. 

We  also  have 


Ir+ 


g{x,y,t)(p{t)dt  = 


dtg{x,y,t)(t){t)  dt 
u*{x,  y,  t)(j){t)  dt 


=  /  u{x,y,t)^{t)  dt 

Jr+ 

where  the  fact  that  dtg{x,y,t)  =  u*{x,y,t)  almost  everywhere  was  used  in  the  second  step 


and  (4.15)  was  used  in  the  third  step.  Thus  we  see  that 

/  {g{x,y,t) -u{x,y,t))^{t)dt  =  0, 


(4.20) 


for  every  0  G  C“(M’'').  Since  g{x,  y,  0)  =  u{x,  y,  0),  we  conclude  from  (4.20)  that  g{x,  y,  t)  = 
u{x,y,t)  almost  everywhere,  completing  the  proof  of  Part  (a). 

Part  (b).  Let  'ip{x,yR)  be  in  x  Y  x  M)^  and  Wperiodic  in  y.  Then  by  using 

integration  by  parts,  we  see  that 


'  Ox]R+ 


dtU^{x,t)  ■ 'Ip  (x, —,t]  dxdt  =  —  u^{x,t)  ■  dt'ip  (x, —,t]  dxdt 

^  JnxR+  \  e  / 

—  j  u^{x,  0)  •  0  -,  oj  dx. 

Sending  e  to  0,  we  obtain 

/  u*{x,y,t)  ■  Ip  {x,y,t)  dxdydt  =  —  u{x,y,t)  ■  dtip  {x,y,t)  dxdydt 

JqxYxR+  JqxYxR+ 


vP’{x,  y)  ■  0  (a:,  ?/,  0)  dxdy. 


>nxY 


(4.21) 
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On  the  other  hand,  using  Part  (a),  we  see  that 


u*{x,y,t)  ■ 'll:  {x,y,t)  dxdydt  = 


'OxyxM+ 


/  dtg{x,  yR)-^)  (x,  y,  t)  dxdydt 

'nxYxK+ 

-  /  g{x,  y,  t)  ■  dt-ip  (x,  y,  t)  dxdydt 
J QxyxR+ 

-  g{x,y,0)  ■'Ip{x,y,0)  dxdy 


>nxY 


>nxYxR+ 


u{x,  y,  t)  ■  dtip  {x,  y,  t)  dxdydt 


u{x,  ?/,  0)  •  '0  {x,  y,  0)  dxdy. 


inxY 


(4.22) 


From  (4.21)  and  (4.22)  we  obtain  that 


>nxY 


{u°{x,  y)  -  u{x,  y,  0))  •  ip  (x,  y,  0)  dxdy  =  0, 


for  every  -0.  Therefore 

u{x,y,0)  =  u'^ix^y), 

almost  everywhere.  Similarly  we  can  show  that 

dtu{x,y,0)  =  n°(x,2/), 

almost  everywhere,  completing  the  proof  of  Part  (b). 


□ 


Proof  of  Lemma^J^  Part  (a).  Since  Al  =  —  ^l,2,  we  will  compute  the  two-scale 

limits  of  Al^iu^  and  then  combine  them  to  show  that  as  £  — >■  0, 

Alu^{x,  t)  ^  Blu{x,  y,  t).  (4.23) 

Let  0  G  C'“(]R^  X  F)^  such  that  0(x,  y)  is  F-periodic  in  y,  and  0  G  C')("(]R’'').  Then  from  the 


dehnition  of  Al,i,  equation  (3.1),  we  see  that 


'  OxM+ 


Al^iU^(x,  t)  ■  tp  —  j  (pit)  dxdt 


f  f  .  (x  —  x)  (x  —  x)  , 

/  /  Xix  —  x) - — - - -  u  ixR)  dx  ■  p)  [x^ —  ]  (pit)  dxdt^ 

InxR+JH^ix)  \x  —  x\'^  \  e/ 


(4.24) 


Since  u^ix,t)  uix,y,t),  we  obtain  using  Proposition  4.9  that,  as  e  — >  0, 


uix,y,t)dy  weakly  in  L^(f2  x  (0, T))" 


'Y 


(4.25) 
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It  follows  from  (4.25)  that,  for  fixed  x, 


lim 


'm+  J  Hj{x) 


Ax  —  xj  ®  (x  —  xj  ,  ,,,  , 

A(a;  —  x) - — - rr - -  u  (x,  t)4>{t)  dxdt 


\x  —  X 


IR+  Jh^{x) 


X{x  —  x) 


[X  —  X) 


(X  —  X) 


\x  —  x\ 


u{x,  y  ,  t)  dy'  )  (j){t)  dxdt 


'Y 


(4.26) 


We  note  that  by  replacing  v{x)  with  u^{x,t)  in  (3.17),  we  obtain 


/  Ax  -  X)  ®  {X  -  X, 

/  X{x  —  X) - - u[x,t)dx 

'  H-y{x) 


\x  —  xP 


^  AjQax 


^  A  max 


'Hs[x) 


'  Hs{x) 


\X  —  X\ 

1 

\x  —  x\ 


dx 


-  dx 


I/s' 


I/s' 


|m^(x,  f)p  dx 


iHsix) 


\u  ||l°o([0^t];L'>(0)3)- 


1/s 


(4.27) 


From  Proposition  3.1,  ||'U^||L°°([o,r];L'>(o)3)  is  bounded.  Thus  from  (4.26),  and  (4.27)  and  by 
using  Lebesgue’s  dominated  convergence  theorem,  we  conclude  that  the  convergence  of  the 
sequence  of  functions  in  (4.26)  is  not  only  point-wise  in  x  convergence  but  also  strong  in 
L^{VtY.  Therefore  we  can  use  Proposition  4.2  and  (4.26)  to  evaluate  the  limit  of  (4.24)  as 
e  — >•  0.  We  hnd  that 


lim 

£—>0 


QxR+ 


t)  ■  -ijj  ix,  -  \  (j){t)  dxdt 


X 


X{x  —  x) 


'  f2xM+  J  H-y{x) 


[X  —  X)  (X  —  X) 


X  —  x\ 


>Y 


u{x,  y',  t)  dy'  )  dx  ■  [x,  —  j  (pi/t)  dxdt, 


X 


Next  we  evaluate  the  two-scale  limit  of  Al^2U^-  We  recall  from  (3.2)  that 
AL,2U%X,t)  =  / 


J  H-y  (X) 

from  which  immediately  follows  that  as  £  — >  0, 


X{x  —  X) - j-x - - dx  U  {x,t), 


X  —  x\ 


Al,2U^ 


'  H Ax) 


Ax-x)®{x-x) 

A(a:  —  x) - j-y — yy - dx  u[x,  y,  t). 


X  —  X 


(4.28) 


(4.29) 


(4.30) 


Combining  equations  (4.28)  and  (4.30),  the  result  (4.23)  follows. 

The  fact  that  the  two  operators  Bl  and  Bs  are  linear  and  bounded  on  the  Banach  space 
Lp^riy))  can  be  shown  by  arguments  similar  to  those  used  in  the  proof  of  Proposition 

EH 

Part  (b).  Since  A^  =  A^^^  —  Ag2,  we  will  compute  the  two-scale  limits  of  A^^u^^  and 
74^2'*^^)  then  combine  them  to  show  that  as  £  — >■  0, 


A%u‘(x,t)  Bsu(x,y,t). 


(4.31) 
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Let  ip{x,yR)  =  where  V’2  e  'tpi  e  and 

0  G  C'“(]R’'').  Then  by  using  (3.13),  replacing  v{x)  with  u'^{x,t),  we  see  that 

t)  ■  'll)  j  dxdt 


'f2x]R+ 


'  £7xR+  J Hs{0) 


X  X 


a  +  Z 


e  e 


z  ^  z 


u^{x  +  ez^  t)dz  ■  Ip  ,  t  j  dxdt. 


(4.32) 


We  recall  that  a  (|,  |  +  z)  is  dehned  by  equation  (1.5).  Without  loss  of  generality,  we  may 
assume  that  a  (|,  |  +  ^)  is  given  by 


X  X 


e  e 


a{-,-  +  z]  =Xi{-]Xf{-  +  z 


x\  fx 


e  J  \e 


Thus  after  a  change  in  the  order  of  integration  in  the  right  hand  side  of  equation  (4.32),  we 
see  that 


'  OxR+ 


'^5(0) 


t)  ■  Ip  (^x,  dxdt 

[  Xf  (-)  Xf  (-  +  t)-z  ipi  f-)  -z  'ip2(x)(p(t)  dxdtdz. 

./ovw+  V £■  /  \£  )  \e  ' 


(4.33) 


Now  we  focus  on  evaluating  the  limit  as  £  0  of  the  inner  integral  in  (4.33).  By  the  change 

of  variables  r  =  a;  +  we  obtain 


Xf  (  -  )  Xf  (  -  +  u^{x  +  ez,t)-z  pJi  (  -  )  'ip2{x)(p{t)  dxdt 

/oxR+  ^ 

Xn{r  -  ez)  xt  Xf  t)-z'ipiQ-  z^ -z  'ip2(r  -  ez)(p(t)  drdt 


(4.34) 


:=  a%z), 

where  xu  denotes  the  indicator  function  of  f2.  We  will  show  that  for  z  G  Hs{0), 


bma%z)  =  /  Xi{y  -  z)xf{y)u{r,y,t)-z  iJiiy  -  z)-z 'ip2{r)(p{t)drdydt. 

JnxYxR+ 

(4.35) 

To  see  this,  we  approximate  xn  by  smooth  functions  Cn  such  that  as  n  — >•  cxo,  Cn(’")  — ^  Xo(?") 
almost  everywhere  and  Cn  — ^  Xn  in  with  1/s  +  l/s'  =  1.  Then  by  adding  and 
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subtracting  —  ez)  to  and  from  xn{r  —  £z)  in  (4.34),  we  obtain  that 


a%z)  =  + 

where, 

:=  /  {xn{r  -  £z)  -  (n(r  -  ez))  X 

Jm.^x'R+ 


n,e  /  \ 

aX  (z)  := 


(e  ~  (D  (e  ~  ~  drdt, 

/  C,n{r-£z)x 
dM3xM+ 

Xf  Xf  t)-z  i>i  (^  -  z^  ■  z  'rp2{r  -  £z)(j)(t)  drdt. 


From  (4.37)  and  by  using  Holder’s  inequality,  we  see  that 


\ar{z)\  < 


\Xn{r  -  ez)  -  Cn{r  -  £z)\'  dr  ]  x 


l/s' 


l/s 


(4.36) 


(4.37) 


(4.38) 


0(t)  dt.  (4.39) 


We  note  that  the  second  term  on  the  right  hand  side  of  (4.39)  is  bounded  above  uniformly  in 


e.  This  follows  from  Holder’s  inequality  applied  to  the  inner  integral  and  the  fact  that  (n^)£>o 
is  bounded  in  L)^^(M+;  L^(fl)^).  On  the  other  hand,  by  the  change  of  variables  r'  =  r  —  ez, 


the  hrst  term  on  the  right  hand  side  of  (4.39)  becomes 


Ixo(’^0  -  Cn(r')r  dr 


l/s' 


which  goes  to  zero  as  n  — >  oo.  From  these  two  facts  and  (4.39),  we  conclude  that  for  all 
£  >  0  and  z  G  Hs{0), 


lim  ai’^{z)  =  0.  (4.40) 

n— >oo 

Now  for  hxed  n,  since  (n  and  'ip2  are  smooth  functions,  we  see  that  as  £  — >•  0,  (C„(r  — 
ez)'4)2{r  —  ez)  — >  Cn(^)t/’2(^)  uniformly.  Therefore,  we  see  from  (4.38)  that 

lim  al’^{z) 


e^O 

=  lim 


R2xR+ 


Cn{r)  Xf  Xf  (^)  U%r,t)-z  ^pl  -  z'^ -z  ^p2{r)(p{t)  drdt 


Cn{r)xi{y  -  z)  Xf  {y)u{r,y,t)-z  tjji  {y  -  z)-z  'ip2{r)(j){t)  drdydt, 


^xVxII 


(4.41) 
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where  in  the  last  step  the  fact  that  (m^)£>o  two-scale  converges  to  u{r,y,t)  was  used.  By- 
taking  the  limit  as  n  — >  cx)  in  (4.41),  we  obtain 


lim  lima2’^(2:) 

n—^oo  £— 

=  /  Xf{y- z)-z'il:2{r)(j){t)drdydt. 
J  flxY  xR+ 


(4.42) 


From  (4.40)  and  (4.42)  and  since 


lima^(^)  =  lim  lim(a”’^(^)  +  a2’^(z)), 

e^O  n— >oo  £—>0 


equation  (4.35)  follows 


From  (4.33)  and  (4.35),  and  by  using  Lebesgue’s  dominated  convergence  theorem,  we 
obtain 


lim 

£— ‘■O 


QxM+ 


Ag  iU^{x,  t)  ■  'tp  “5  dxdt 


=  r-f^  /  Xf{y  -  z)Xi{y)u{r,y,t)-z -ipiiy  -  z)-z 'Ip2{r)4>{t)drdydtdz 

JHsiO)  JQxYxm+ 

=  I  I  /  Xi{y- z)xi{y)u{r,yR)-zi)i{y- z)-zdydzi)2{r)(t){t)drdt, 

JnxR+  Jhs(o)  \^\  Jy 

(4.43) 

where  we  have  changed  the  order  of  integration  in  the  last  step.  After  shifting  the  domain 


of  integration  in  the  inner  integral  of  the  right  hand  side  of  equation  (4.43),  we  obtain 


Xf  {y  -  z)  Xf  (y)  u(r,  y,t)-z  pJi  {y-  z)-zdy 


'Y 


Xi  {y)  Xf{y  +  z)u{r,y  +  z,t)-z  tpi  {y)-zdy 


lY-z 


=  y  Xf  (y)  Xf  {y  +  z)  M(r,  y  +  z,  t)-z  {y)-zdy, 


(4.44) 


where  in  the  last  step  the  fact  that  the  integrand  is  F-periodic  in  y  was  used.  Substituting 


(4.44)  in  equation  (4.43),  then  by  changing  the  order  of  integration  we  obtain 


lim 

£— ‘■O 


QxR+ 


A^g  -pf{x^  t)  ■  Ip  ^x,  — ,  t  j  dxdt 


z  z 


inxR+  Jy  J Hs{o) 

/  / 

'nxYxm+  JHsiy) 


Xf  (y)  Xf  {y  +  z)  -—^u{r,  y  +  z,  t)dz  ■  ipi{y)dy  ip2{r)(p{t)drdt 

pI 

Xi{y)Xi  {y)  — — — ~  “(a  drdydt. 


(4.45) 


In  the  last  equality  the  change  of  variables  y  =  y  +  z  was  used. 


Bacim  Alali  and  Robert  Lipton 


32 


Next  we  evaluate  the  two-scale  limit  of  Ag2'^^-  Let  'ijj  he  a  test  function  in  J'.  Then  by 
using  (3.16),  replacing  v{x)  with  u^{x,t),  we  obtain 


'  f2x]R+ 


f)  •  '?/’  ,  t  j  dxdt 


'  OxE+  J  Hs{0) 


X  X 


a  +  z 


e  e 


z  ^  z 


dz  u^{x^  t)  ■  '0  dxdt. 


The  right  hand  side  of  (4.46),  after  changing  the  order  of  integration,  is  equal  to 

I  ^  [  «  f— , — \- z^  u^(x,t)-z  (x, —,t] -z  dxdtdz. 

JHsio)\^r  JnxR+  e  J  \  e  / 

Using  the  fact  that  (M^)e>o  two-scale  converges  to  u{x,y,t),  we  see  that  for 

2  e  HsiO), 


lim 

£—>0 


T2xR+ 


a  ,  — \-  z^  u^{x,  t)-z  -z  dxdt 

a  {y,y  +  z)  u(x,  y,t)-z  (x,  y,  t)-z  dxdydt. 


'nxYxR+ 


(4.46) 


(4.47) 


(4.48) 


From  (4.46),  (4.47)  and  (4.48),  and  by  using  Lebesgue’s  dominated  convergence  theorem,  we 
obtain 


lim 


T2xR+ 


Ag  2'^^{x,  t)  ■  'ifj  (^x,  — ,  t  j  dxdt 


iHsiO) 


—r^  /  a{y,y  +  z)u{x,y,t)-z  Ip  {x,yA)-zdxdydtdz 


lnxYxR+ 


(4.49) 


By  changing  the  order  of  integration  and  then  using  the  change  of  variables 
y  =  y  +  z,  we  conclude  that 


lim 

£— >0 


f^xR+ 


Ag  2'^^{x,  t)  ■  Ip  — ,  f  j  dxdt 


'nxYxR+  JHsiy) 


a  {y,  y)  — — — —  dy  u{x,  y,t)  -  ip  {x,  y,  t)  dxdydt. 


(4.50) 


Equation  (4.31)  follows  from  combining  (4.45)  and  (4.50),  completing  the  proof.  □ 


5  The  Macroscopic  Equation  and  Downscaling 

The  aim  of  this  section  is  to  justify  the  main  results  of  Section 
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5.1  Derivation  of  the  Macroscopic  Equation 

We  begin  this  section  with  the  following  observation.  Let  0  be  a  function  in  Then 


a{y,yA  Wv)  -  ■#>(!/))  =  0. 


(5.1) 


lYJHsiy)  \{y~y)\^ 

To  see  this,  we  note  that  using  Fubini’s  theorem  and  the  assumption  that  0  is  F-periodic, 


the  double  integral  in  (5.1)  can  be  written  as 


V  JHsiy) 


\{y-yW 


Y  JHsiy)  \{y-y)r 


(5.2) 


where  in  the  last  equality  we  have  used  the  fact  a{y,y)  =  a{y,y).  Comparing  the  double 


integral  in  (5.1)  with  (5.2)  the  result  follows. 

Now  let 

u^{x,t)  =  j  u{x,y,t)dy. 

Then  from  Proposition  4.9,  we  have  that  u^{x,t)  is  the  weak  limit  of  u^{x,t)  in  LP{fl  x 
(0,T))^.  To  identify  the  equation  that  solves,  we  integrate  (4.11)  over  Y  to  obtain 


d‘^u^{x,t)  =  /  X{x  —  x) 

J  Hjix) 


[X  —  xj  <Y>  [X  —  xj 
\x  —  x\^ 


{u^{x,  t)  —  {x,  t))  dx 


Hi 


Y  JHsiy) 


a{y,  y)  — — d,  t)  “  y,  t))  dydy  (5.3) 


+  j  b{^,y,t)dy. 


Using  (5.1),  the  second  integral  on  right  hand  side  of  (5.3)  is  equal  to  zero  for  all  x  G  fl  and 
t  E  (0,T).  Thus  solves 


d^u^{x,t)=  /  X{x  —  x)- 


JH-yix) 

supplemented  with  initial  data 


\x  —  x\ 


{u^  {x,t)  —  {x,t))  dx  +  J  b{x,y,t)dy,  (5.4) 


u  (a;,0)  =  J  u^{x,y)dy,  dtu  (x,  0)  =  j  v^{x,y)dy. 


(5.5) 


The  initial  value  problem  (5.4)-(5.5)  can  be  written  as  the  following  operator  equation  in 

UiSlf 

u^{t)  =  ALU^{t) +  h{t),  te[0,T] 

M'^(O)  =  -u°, 

'U'^(O)  = 


(5.6) 
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b{x,t)  = 

/  K^,y,t)dy, 

Jy 

u°(x)  = 

/  u^{x,y)dy,  and 

Jy 

v°(x)  = 

j  v°{x,y)  dy. 

We  have  seen  from  the  proof  of  Proposition  |3.1  that  Al  is  a  bounded  linear  operator  on 


thus  G  C‘^{[b,T]]  U{VtY)  is  the  unique  solution  of  5.6 


To  complete  the  proof  of  Theorems  2.1,  2.3,  2.5,  and  2.7,  we  show  that  is  in 
C^([0,T];  (7(12)^),  when  the  initial  data  and  are  in  C'(r2)^,  and  the  loading  force  b  is 
in  C([0,T];  (7(12)^).  In  fact,  it  suffices  to  show  that  the  linear  operator  Al  is  bounded  on 
the  Banach  space  of  continuous  functions  equipped  with  the  uniform  norm.  So  we  let 

V  G  C(f2)^  and  denote  the  uniform  norm  on  C(f2)^  by  ||  ■  || (7(0)3-  Then,  we  recall  from  (3.5) 
that  Al  =  Al^i  +  ^l,2,  where  Al^i  and  Al^  can  be  written  as 


Al,iv{x)  =  [  v{x  +  Od^, 

Jh^{o)  Is  I 

Kl 


(5.7) 

(5.8) 


respectively.  Taking  the  norm  in  (5.7)  we  see  that 


Pn,w|lc(o)3  =  max 


< 


< 


max  A(7)  max  /  —  |n(a:  +  7)|(i7 

€er/,{0)  xeo  7^,(0)  1^1 


inax  A(0  /  ll'y||c(o)3- 

€er/^(o)  J  Jh^{q)  Is  I 


Thus  Al,i  is  bounded  on  C(f2)^.  It  is  clear  that  Al,2  is  also  bounded  on  (7(12)^,  and  therefore 
Al  is  bounded  completing  the  argument. 


5.2  Justifying  the  Downscaling  Step 

In  this  section  we  prove  Theorems  BIHIll  and  |2.8[  We  begin  by  showing  that  for  fixed 
t  G  (0,T), 
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By  shifting  the  domains  of  integration,  equation  (4.11)  can  be  written  as  follows 

/ 

iH-yiO)  ISI 


d^u{x,  y,  t)  = 


uix  +  C,  y\  t)  dy'  -  u{x,  y,  t) 


Y 


z  z 


+  /  a(2/,  y  +  z)—^  {u{x,  y  +  zR)-  u{x,  y,  t))  dz 

Jhs{0)  Bh 

+  b{x,yR). 


(5.9) 


Since  u{x,  y,  t)  is  in  Q  and  solves  (5.9)  with  initial  conditions  (4.12)  and  (4.13),  then  u{x,  -,  t) 
is  in  C^([0,T];  and  solves 


S,\(i,  pi)  = 


/ 

'H^iO)  Isl 


'Hs{0) 


XX  \  Z  Z 

“h'l  +  Vl^F 


u{x  +  ^,y,t)dy-u(^x,^,t^'^  d^ 

(^u(^x, — \- z,t^  —  u  (^x, dz  (5.10) 


+  b(^x,  pt)  , 

supplemented  with  initial  conditions 

u(x,y,0)  =  tl“(a:,  j), 
d,u(x,y,0)  =  li“(a:,  j). 


(5.11) 

(5.12) 


We  let  e^{x,t)  =  u^{xR)  —  u{x,  -,t).  Then  by  subtracting  (5.10)  from  (1.10),  we  hnd  that 
G  C^([0,T];  L®(fl)^)  solves 


d‘le^{x,t)  =  A^e^{x,t)  +  d^{x,t), 
e^(a;,0)  =  0, 
dte^{x^b)  =  0. 


where  A^  is  given  by  (3.7)  and  d‘^{x,t)  is  given  by 


d^{x,t)  =  d'i{x,t)  +  dg{x,t), 


(5.13) 

(5.14) 

(5.15) 


(5.16) 


dl{x,t)  =  /  A(0 


lep 


U  (  x  + 


u{x  +  ^,y'A)dy' 


>Y 


(5.17) 


d‘j(x,t)  = 


/  /XX  \  Z  z 

'hs(o)  ^  \z\ 


(^u  (^x  +  ez,  —  +  z,t^  —  u  ^a;,  —  +  z,t^'^  dz. 


(5.18) 


Since  A^  is  bounded,  the  solution  of  (5.13)-(5.15)  is  explicitly  given  by 


pt  ^ 

eAxR)  =  I  Y1 


(t-r) 


2n+l 


r2n  +  Id 


(7h)"d"(a:,r)  dr. 
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Thus 


(f  _  ^\2n+l 


n=0 


< 


Vm 


sinh  (yM {t  -  t)^  \\d‘'{-,r)\\Ls^Qp  dr  (5.19) 


where  in  the  second  inequality  we  have  used  the  fact  that  is  bounded  above  by  an  M  >  0 
independent  of  e. 

In  the  following  sections  we  will  show  that  for  t  G  (0,T), 


lim  ||ci^(-,f)||i.(Q)3  =  0, 


(5.20) 


for  each  of  the  four  cases  of  initial  and  loading  conditions  that  has  been  introduced  in  Section 
On  the  other  hand,  from  (5.16)-(5.18)  and  the  fact  that  u  is  continuous  on  [0,  T],  it  follows 
that  d^{-,T)  is  continuous  on  [0,f]  for  t  <T.  Thus,  from  equations  (5.19)  and  (5.20),  and 
Lebesgue’s  convergence  theorem,  we  see  that 


lim  ||e''(-,f)||Lqn)3  =  0, 


from  which  the  result  follows. 


In  order  to  prove  (5.20),  we  will  make  use  of  the  following  observation: 

The  solution  of  each  cell-problem  of  Section  has  zero  average  over  the  unit  cell.  To  see 


this,  we  integrate  equation  (|2.5|)  over  Y  to  obtain 
f{t)  = 


supplemented  with  initial  conditions 


dydy-K  r(t),  (5.21) 


r(0)  =  0,  r(0)  =  0. 


(5.22) 


Using  (5.1),  the  integral  on  the  right  hand  side  of  (5.21)  is  equal  to  zero  for  all  t  G  (0,T). 
Thus  f  solves 

f{t)  =  —K  f{t),  (5.23) 


supplemented  with  zero  initial  conditions.  Obviously  the  solution  of  (5.23)  is  given  by 


'■(</.*)  %  =  »■(*)  =  0. 


(5.24) 


lY 


for  all  t  G  (0,T).  Similarly  we  can  show  that 


r^{y,t)dy  =  Y{t)  =  0, 


(5.25) 


'Y 


for  all  t  G  (0,T),  where  Y  is  the  solution  of  (2.12),  (2.19),  or  (2.26). 
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5.2.1  First  Case 


In  this  section  we  complete  the  proof  of  Theorem  2.2  by  showing  that  equation  (5.20)  holds 
true  when  b,  vP,  and  are  given  by  (2.2).  We  also  prove  the  error  estimate  (2.8). 

Using  the  fact  that  r{y,t),  the  solution  of  the  cell  problem  (2.5)-(2.6),  has  zero  average 
over  Y,  and  by  linearity,  it  is  easy  to  check  that  u^{x,t)  +r{y,t)  solves  (4.11  )-(4.13 ),  where 
is  the  solution  of  (2.3)-(2.4).  Thus  by  uniqueness  we  conclude  that 


uix.yR)  =  u^{xR)  +r{yR). 


(5.26) 


Using  this  representation  of  u{x,y,t)  and  from  equations  (|5.17)  and  (5.18),  we  see  that 
d^^{x,t)  and  dg{x,t)  are  now  given  by 


dl{x,t)  = 


.  fx  +  p 


/  /XX 

/  a  +  z 

Ihs{0)  £ 


e 

z  ^  z 


,tj  de,  (5.27) 

(u^ {x  +  ez,t)  —  {x,t)^  dz,  (5.28) 


respectively. 


Changing  variables  of  integration,  equation  (|5.27|)  becomes 
dl{x,t)  =  [ 


'  [x) 


Ax-x)®{x-x)  fx 
X{x  —  X) - - rl  -,f)  dx. 


\x  —  xP 


(5.29) 


Since  r{y,t)  is  U-periodic  in  y  and  from  Proposition  4.10,  we  see  that  for  fixed  f,  as  £  — >  0 

^  /  r{yA)  dt  =  b  weakly  in  L\VLf. 


'Y 


Thus  from  (5.29)  we  obtain  that 

limd|^(x,  t)  =  0, 

for  x  G  12  and  t  G  (0,T).  It  follows  from  Lebesgue’s  convergence  theorem  that 

lim||di(-,t)||Lqo)3  =  0, 

£—►0 


(5.30) 


for  t  G  (0,T).  On  the  other  hand,  by  taking  the  Euclidean  norm  of  dg{x,t)  in  (5.28),  we 
obtain 


|d|(a;,f)|  < 


'Hs{0) 


—  \u^{x  +  ez,  t)  —  u^{x,  t)  I  dz, 


(5.31) 


where  Omax  =  max  a{y,y').  Since  G  C'^([0,T];  C(12)^)  (see  Section 

y,y'& 

for  X  G  12  and  t  G  (0,  T) 


5.1),  it  follows  that 


lim  |d|(x,f)|  =  0.  (5.32) 
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Thus  using  Lebesgue’s  convergence  theorem,  we  obtain 


lim||4(-,t)||i.(f^)3  =  0, 


(5.33) 


for  t  e  (0,T).  Equation  (5.20)  follows  from  equations  (5.30)  and  (5.33) 


Now  we  prove  the  error  estimate  (2.8).  By  setting  A  =  0  in  equation  (2.3),  we  see  that 
its  solution  is  given  explicitly  by 


(xR)  =  uo{x)  +  t  vo{x)  +  I  {t  —  t)1{x,t)  dr. 


(5.34) 


By  assumption  uq,  Uq,  and  are  in  C°’^(r2).  Thus  for  2;  e  ^^5(0),  we  see  from  (5.34)  that 


\u^ {x  +  ez,t)  —  {x,t)\  <  C\ez\^  +  t  C\ezp  +  /  (t  —  T)C\ezp  dr 

Jo 

=  C  (l  +  t+^-]\z\^e>^,  (5.35) 


for  some  C  >  0.  We  use  this  bound  in  inequality  (|5.31)  to  obtain 

f  \z\^-Uze^. 


\dsixR)\  <  C  1  +  t  +  -  a 


^‘2 


>Hs{0) 


(5.36) 


Since  A  =  0  we  see  from  (5.16)-(5.18)  that  d^  =  dg.  Therefore  from  (5.36),  after  a  simple 
calculation,  we  obtain 


^/3+2  /  4.2 ' 

IM'(-.*)IIl-(!1)s  <  (  1  +  (  +  -  )  ef 


(5.37) 


By  using  (5.37)  to  bound  ||(i^(-,  r)||£,s(Q)3  in  (5.19),  the  error  estimate  (2.8)  follows. 


5.2.2  Second  Case 


In  this  section  we  complete  the  proof  of  Theorem  2.4  by  showing  that  equation  (5.20)  holds 


true  when  b,  vP,  and  are  given  by  (2.9).  We  also  prove  the  error  estimate  (2.15). 


Using  the  fact  that  r(j/,t),  the  solution  of  the  cell  problem  (2.12)-(2.13),  has  zero  average 


over  U,  and  by  linearity,  it  is  easy  to  check  that  u^{x,t)  +  {yR)hj{x)  solves  (4.11)- 


(4.13),  where  is  the  solution  of  (2.10)-(2.11 ).  Thus  by  uniqueness  we  conclude  that 


uix,y,t)  =  u^{x,t)  +  ^r\y,t)hj{x). 

i=i 


(5.38) 


Using  this  representation  of  u{x,y,t)  and  from  equations  (5.17)  and  (5.18),  we  see  that 
d\^{xR)  is  now  given  by 


dl{xR) 


i=i 


x  +  ^ 


,t)  hj{x  +  ^)d^, 


(5.39) 
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and  dg{x,t)  can  be  written  as 


ci|(x,  t)  =  t)  +  t), 


(5.40) 


where, 


d^s,i{x,t)  = 

d%,2{xR)  = 


'Hs(Q) 


/XX  \  Z  ®  Z  ,  TT  ,  ,  W/N\,  /.N 

a  (  — \-  z]  iu^ix  +  ez^  t)  —  (x,  t))  dz,  (5-41) 

\e  e  /  Izr 


/  /XX  \  z  z 
/  a[-,-  +  z] 

iHsiO)  ^  ’  kh 


^ — \-  z,t^  {hj{x  +  ez)  —  hj{x))  dz. 
i=i 


(5.42) 


Applying  the  methods  developed  in  Section  5.2.1  for  (5.30)  and  (5.33),  we  can  show  that 
for  t  e  (0,  T), 

lim  ||ci|(-,f)||i.(n)3  =  0,  (5.43) 


and 


lim||d|,i(-,t)||L.(o)3  =  0. 


It  remains  to  show  that  for  t  G  (0,  T), 


hm||42(-)^)iUdf^)3  =  0- 
£— >0  ’ 


(5.44) 

(5.45) 


From  eqnation  (5.42),  we  see  that 


^  “max  /  +  \hj{x  +  ez)  -  hj{x)\dz,  (5.46) 

I  u.rn\  IZl  I  \E  / 


i=i 


where  Omax  =  max  aiy^y').  Since  |  <  s  <  cx),  we  can  choose  s',  with  |  <  s'  <  cx),  and  s". 


y,y'eY 


with  1  <  s"  <  3,  snch  that  1/s  +  1/s'  +  1/s"  =  1.  By  Holder’s  ineqnality  we  obtain 


\ds,2ix,t)\  <  ar 


Hs{0) 


1/s"  3 

dz) 

j=i 


Hs{0) 
1/s 


X 


r  — V  zR 


dz 


I/s' 


x(  I  \hj{x  +  ez)  —  hj{x)\^  dz 

Hs{0) 


It  is  easy  to  see  that 


Hs{0) 


;  /X 


— \-  z,t 


l/s' 


dz]  <  ||r^(-,t)||LV(o)3. 


(5.47) 


(5.48) 


Thns  from  (5.47)  and  (5.48),  and  by  using  the  triangle  inequality  in  L^,  we  obtain 

l/s"  3 


II4,2(-,^)II  <  “r, 


Hs(0) 


X 


dz]  ||r^(-,t)||LV(f^)3 

i=i 

l/s 


I hj {x  +  ez)  —  hj (x)  dzdx 
Q  JHs(0)  J 


\  i/« 


(5.49) 
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Since  hj  is  continuous  on  we  obtain  from  Lebesgue’s  convergence  theorem  that 


lim 

£•— »'0 


O  JHs{0) 


\hj{x  +  ez)  —  hj{x)\^  dzdx  =  0. 


(5.50) 


Equation  (5.45)  follows  from  (5.49)  and  (5.50).  This  shows  that  (5.20)  holds  true  for  this 
case. 


Now  we  prove  the  error  estimate  (2.15).  By  setting  A  =  0  in  equation  (2.10),  we  see  that 
its  solution  is  given  explicitly  by 


u^{xR)=  /  {t  —  t)F{t)  dr  h{x). 


(5.51) 


By  assumption  h  is  in  Thus  for  G  ^^^(O),  we  see  from  (5.51)  that 

rt 


\u  {x  +  ezR)—u  {xR)\  <  C\ez\^  I  {t  —  T)F{T)dT, 

Jo 


(5.52) 


(5.52),  we  see  that 


for  some  C  >  0.  Taking  the  Euclidean  norm  in  both  sides  of  (5.41)  and  using  the  bound 

(5.53) 


|d|i(x,f)|  <  /  (f-r)F(r)dr  /  ^dze^, 

Jo  JHs(0) 


and  it  follows  that 


^  {t-T)F{T)dT^  .  (5.54) 


On  the  other  hand  from  (5.49),  after  a  straight  forward  calculation,  we  obtain 

l/s' 


^  C'ctmax  dvT 


53- 


3 -s' 


A^/3+3  \  JL 

47r|0|  ^  )  5^||r’^(-,t)||L«'(0)3^^- 

i=i 


s/3  +  3 


(5.55) 


Since  A  =  0  we  see  that  d^  =  d^-^F  d^g2-  Therefore  by  combining  (5.54)  and  (5.55)  to  bound 
\\d^i;r)\\L  qn)3  in  (|5.19|),  the  error  estimate  (|2.15|)  follows. 


5.2.3  Third  and  Fourth  Cases 


Arguments  similar  to  those  presented  in  Section  5.2.2  show  that  equation  (5.20)  holds  true 


when  the  loading  and  initial  conditions  are  given  by  (2.16)  or  (2.23).  Also,  the  proofs  of  the 


error  estimates  (2.22)  and  (2.29)  are  similar  to  the  proof  of  (2.15)  provided  in  Section  5.2.2 


For  completeness,  we  explicitly  provide  the  functions  M2,(t)  and  M^it)  of  Theorems  2.6  and 


2.6,  respectively.  The  function  M2,(t)  is  given  by 


M3{t)  =  J  sinh  (^^/M {t  -  fsir)  dr 
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where 


/  _  XI3+2  /  r3 


X  V-" 


47r 


^s/3+3  \  1/^  3 


s/3  +  3 


i=i 


L='(0.f 


and  +  solves  (2.19)-(2.20). 


The  function  M^if)  is  given  by 


M^{t)  =  /  _ sinh  ( \fM {t  —  r))  f^ir)  dr 

Jo  yM  ^  ' 


where 


/  _  A'^+2  /  \  /  a^/9+3  \  3/®  3 

/4(i)  =  Ca„„|SJ|‘/<|^4^|F|_i+(^4x— j  {*^713^3)  Ell’''(V)ll 


L^'inp 


and  +  solves  ( 2.26 )-( 2.27). 


This  completes  the  proofs  of  Theorems  2.6  and  2. 


6  Fluctuating  Long-Range  Bond  Model 

In  this  section,  we  present  a  new  multiscale  analysis  method  for  computing  the  deformation 
of  hber-reinforced  composites  modeled  by  the  peridynamic  formulation.  This  is  done  for 


the  Fluctuating  Long-Range  Bond  model  described  in  Section  |1.2[  The  method  provides  a 
computationally  inexpensive  multiscale  numerical  method.  This  is  described  by  Theorem 
6.1[  A  homogenization  result  for  this  model  is  expressed  in  Theorem  6.2 


We  begin  by  recalling  the  peridynamic  equation  of  motion  for  this  model.  By  expanding 


in  equation  (1.15),  then  collecting  the  xf  terms,  we  obtain 


d^u%x,t)  =  xl{x)  [  {Cf-eCra)— — — —{u%x,t)-u%x,t))dk 


X  —  X 


(x-x)(^(x-x) 

eCm - - tt; - [u  {x,t)—U  [x,t))dx, 


'  Hs{x) 


\x  —  x\ 


where  the  hrst  integral  in  (6.1)  is  a  line  integral  over  the  set 


Is{x)  =  {x  E  Hs{x)  such  that  x  —  x  is  parallel  to  n}. 

The  initial  conditions  supplementing  this  equation  are  given  by 

u^{x,0)  =  (x) , 

dtu^{x,0)  =  (x) . 


(6.1) 


(6.2) 

(6.3) 


The  well-posedness  of  equation  (6.1  )-(|0)  is  provided  in  Section [6T|  (Proposition [6^ . 
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Theorem  6.1  (Downscaling).  Let  G  C^([0,T];  be  the  solution 

,3  =  0, 


of  (6.1)- (6.3),  where  1  <p  <  oo.  Then  for  t  e  [0,T], 

lim  t)  —  (xAx)w{x,  t)  +  u^{x)  +  tv^{x)) 

£—>0  "  '  ' 

where  w  G  C^([0,T];  Lp{LI)^)  is  the  solution  of 
d(w{x,t)  =  I 


\LP{np 


(6.4) 


[X  -  X  {X  -  X 

Cf - - tt; - {w{x,  t)  —  w{x,  t))  dli 


\X  —  X\ 


+ 


{x-x)®{x-x)  ,  ^.,0/ 


lq{x) 


\x  —  x\ 


{w{x)  +  tv^{x)  —  {w{x)  +  tv^{x)))  dli 


supplemented  with  the  initial  conditions 


w{x,b)  =  0, 
9w(a;,0)  =  0. 


Moreover,  fort  G  [0,T]  the  error  in  (6.4)  is  estimated  by 


I'x’ix.t)  -  {x)(x)w{x,t)  +!1°(I')  <  e 


(6.5) 


(6.6) 

(6.7) 


(6.8) 


where 


M^{t)  =  (  \\u°\\LP{n)3Cosh^/Mt  +  ||f°||LP(n)3-^  sinh  a/m^ 


and  where  M  is  a  positive  constant. 


Theorem  6.1  is  proved  in  Section  6.2 


The  macroscopic  peridynamic  equation  for  this  model  is  given  by 


d(u^{x,t)  =  /  C{ 
Jq{x) 

+  {Of  - 1)  /  a 

Jl«(x) 


.X  —  X)  (X  —  X 


\x  —  x^ 

[x  —  x)  ®  [x  —  X 


\x  —  X 


2  yU^{x,t) -U^{x,t))dh 

U^{x)  +  tv^{x)  —  (u^(x)  +  tv'^(x)))  dlx, 


supplemented  with  initial  conditions 


u^{x,0)  =  (x) , 

dtu^{x,0)  =  (x) . 


(6.9) 


(6.10) 

(6.11) 


Here  the  macroscopic  displacement  is  the  weak  limit  of  the  sequence  of  displacements  u^. 
This  is  described  by  the  following  theorem. 
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Theorem  6.2  (Homogenization).  Let  E  C^([0,T];  be  the  solution 


of  (6.1)- (6.3),  where  I  <  p  <  oo.  Then  for  t  E  [0,T],  as  e  ^  0, 

M^(-,  t)  — >  t)  weakly  in 


where  E  C^([0,T];  L^ikl)'^)  is  the  solution  of  (6.9)-(6.11).  Equivalently,  can  be  com¬ 
puted  as  follows 

u^{x,  t)  =  9fw{x,  t)  +  u^{x)  +  tv^{x),  (6.12) 


where  w  solves  (6.5)-(6.1). 


Theorem  6.2  is  proved  in  Section  6.2 


Remark  6.  We  observe  that  the  macroscopic  peridynamic  equation  (6.9)  has  a  nonzero 


loading  force,  although  the  original  peridynamic  equation  (6.1)  has  no  loading  force.  The 


physical  interpretation  for  this  phenomenon  is  not  well-understood  up  to  this  point. 

6.1  Existence  and  Uniqueness  Results 

Without  loss  of  generality,  we  may  choose  the  hber  direction  to  be  parallel  to  the  a:i-axis.  So 
let  n  =  (1,  0, 0).  We  note  that  the  matrix  multiplying  (pf{x,  f)  —  u^{x,  t))  in  the  hrst  integral 


of  (6.1)  is  now  given  by 


(X  —  X)  [X  —  X 


\x  —  x\ 


for  xi  7^  Xi.  Thus  equation  (6.1 ),  after  shifting  the  domain  of  integration  in  the  hrst  integral, 
becomes 

rS 


d^u^xR)  =  {Cf  -  eC,^)  x({x)  J  {ul{x -E  (1,0,0),  t)  -  ul{x,t))  dl 

f  ^  (x  -  x)  ®  (x  -  x)  .  _ 

-|-  /  eCm  - — - rr - (u(x,t)—u(x,t))dx. 

JHsix)  \X  -  X\^ 


(6.13) 


Let  V  =  (vi,V2,V3)  E  LP(Q)^  with  1  <  p  <  cx).  Then  we  dehne  the  following  operators 


A{v(x)  =  Cf  J  {vi{x -E  (1,0,0))  —  vi{x))  dl, 

A^fv(x)  =  Xi(x)  Aiv(x), 

.  ,  ,  f  ^  (x  -  x)  ®  (x  -  x'^ 

A^V(X)  =  /  - -2 2^ - 

Jhs(x)  \x  -  X\^ 

a 


v(x)  —  v(x))  dx. 


=  A(  +  e{A^--^AI 


(6.14) 

(6.15) 

(6.16) 

(6.17) 


The  initial  value  problem  (6.1)-(6.3)  can  be  written  as  the  following  operator  equation  in 

LP(n)^ 


u%t)  =  A^u%t),  t  e[0,T] 

M^(0)  =  u°, 

M=(0)  =  v°. 


(6.18) 


Existence  and  uniqueness  of  solution  of  (6.18)  is  given  by  the  following  proposition. 
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Proposition  6.4.  Let  1  <  p  <  oo.  Then 
(a)  The  operator  is  linear  and  uniformly  bounded  on 


(b)  Equation  (6.18)  has  a  unique  classical  solution  G  ^^([OjT];  which  is  given 

by 


°°  4.2n  °°  j-2n+l 


(v4")V 


^  (2n  +  1)! 

n=0  ^  ^ 


(6.19) 


Proof.  Part  (a).  First,  we  show  that  the  linear  operator  Am  is  bounded  on  Let 

V  G  LP{Q)^.  Then  from  (6.16),  Am  can  be  written  as 


C*m(Am,l  Am, 2)) 


where 


Am,l^^(a:)  = 

Am,2^^(a:)  = 


/  [X  —  X)  [X  —  X) 

Ihs(x)  \x-x\^ 

r  {x  —  x)  ®  {x  —  x) 
Ihs(x)  |h-a;|2 


v{x)  dx, 
dxv(x). 


From  equation  (6.20)  we  see  that 


(6.20) 

(6.21) 


IIAmp^l 


LP(0)3 


< 


|n(a;)|  dx  )  dx 


Jo,  \JH-y(x) 

—  1^1  lkllLP(0)3) 


(6.22) 


where  the  fact  that  ||n||Li(o)3  <  ||'y||LP(o)3  was  used  in  the  last  step.  This  shows  that  Am,i  is 
bounded  on  Lp{Q)^.  The  boundedness  of  Am, 2  is  clear.  Therefore  Am  is  bounded  on 

given 


Next  we  note  that  Af  is  bounded  on  which  is  a  consequence  of  Lemma  6.5 


at  the  end  of  this  section.  Thus  it  follows  from  (6.15)  that  Af  is  uniformly  bounded  on 

LP{n)P 


Combining  these  results  with  equation  (6.17),  it  follows  that  A^  is  uniformly  bounded  on 
L^(f2)^,  completing  the  proof  of  Part  (a). 


The  proof  of  Part  (b)  is  similar  to  the  proof  of  Part  (b)  of  Proposition  3.1 
Lemma  6.5.  Let  v  be  in  L^{Q)^,  where  1  <p  <  00,  and  define 


v{x)  =  J  v{x  +  (1,0,0))  dl. 


□ 


Then  v  is  in  L^(f2)^  and 


h||LP(o)3  <  27  ||n||LP(o)3. 


(6.23) 
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Proof.  From  the  definition  of  v  it  is  easy  to  see  that 

/  \v{x)\^  dx  <  {  \v{xi  +  I,  X2,  Xs)\  dl  \  dxidx2dxs. 


n  \J-5 


(6.24) 


Using  Holder’s  inequality  in  the  inner  integral  with  v  G  and  1  G  where 

1/p  +  1/p'  =  1,  we  obtain 


/  \v{x)\Pdx  <  i26y/P'  /  /  \v{xi  +  l  ,  X2,  Xs)  1^  dl  dxidx2dx^ 
In  JnJ-s 

=  {26YIp'  [  [  \v{xi  +  l  ,  X2,  X3)  1^  dxidx2dx^dl, 

J-5  Jn 


(6.25) 


by  Fubini’s  theorem.  We  extend  v  to  by  setting  v  =  0  outside  H.  Then  by  the  change  of 


variables  xi  =  xi  +  I  in  the  inner  integral  of  (6.25),  we  obtain 


/  \v{xi  +  l,X2,Xs)\P  dxi  <  /  \v{xi,X2,X3)\^  dxi. 

In  Jn 


Using  this  estimate  in  (6.25),  we  conclude  that 

[  \v{x)\Pdx  <  {25ylP'{25)  f  \v{x)\Pdx, 


and  (6.23)  follows,  completing  the  proof. 


(6.26) 

□ 


6.2  Multiscale  Analysis  Using  the  Semigroups  Approach 


The  aim  of  this  section  is  to  prove  Theorems  6.1|  and  |6.2[  Our  approach  is  summarized  by 
the  following  steps: 

1.  Compute  the  two-scale  limit  u{x,  y,t)  of  the  sequence  using  the  explicit  represen- 

(6.27) 


tation  of  u^,  equation  (6.19).  We  show  that  for  fixed  t  G  [0,T],  as  £  — 0, 


u^xR)  u{x,yy), 


where  u  is  given  by 


X2n 

u{x,  y,  t)  =  u^{x)  +  tv^{x)  +  Xf(p)  (Hf)  V(a;) 

n=l  ' 


°°  j.2n+l 

n=l  ^  ' 


(6.28) 


2.  Compute  d^u  in  (6.28)  then  use  it  to  identify  the  two-scale  limit  equation.  We  find 


that  u  G  C^([0,T];  LP{Qy)  uniquely  solves 

dtuix.yR)  =  Aiu{x,yy) +  b{x,y,t), 
u{x,y,0)  =  u^{x), 

du{x,y,0)  =  n°(a;). 


(6.29) 
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where  b  is  given  by 

Kx,  y,  t)  =  ixiiy)  -  +  tv^){x). 

Here  the  operator  Af  is  dehned  as  follows.  For  v  G  Lp{Q  x 


Aiv{x,y)  =  Cf  j  {vi{x  + (1,0,0),  y) -vi{x,y))  dl. 


(6.30) 


3.  The  macroscopic  equation  is  found  by  integrating  (6.29)  over  Y.  We  hnd  that  the 
macroscopic  displacement  solves 


dfu^(x,t)  =  AfU^  (x,t) +  b(x,t), 
u^(x,0)  =  u^(x), 

du^(x,0)  =  v^(x), 

where  b  is  given  by 

b(x,  t)  =  (Of  —  l)y4f(M°  +  tv^)(x). 
Here  for  hxed  f  G  [0,  T],  as  £  — >  0, 

u%-,t)  ^u^(-,t)  weakly  in  LP(n)^. 


(6.31) 


(6.32) 


4.  The  two-scale  limit  u  can  also  be  computed  by  the  following  method.  This  method  is 
numerically  inexpensive. 


u(x,  y,  t)  =  Xi(y)w(x,  t)  -1-  M°(a;)  tv^(x), 

where  w  G  C^([0,T];  Lp(Vl)'^)  solves 

dfw(x,t)  =  A{w(x,t)  +  A{(u^  +  tv^)(x), 
w(x,0)  =  0, 

dw(x,0)  =  0. 


(6.33) 


(6.34) 


It  follows  from  integrating  (6.33)  over  Y  that  can  also  be  computed  by 


.H, 


u 


X,  t)  =  9fw(x,  t)  -I-  u^(x)  +  tv^(x). 


(6.35) 


5.  Extend  u  by  periodicity  from  hi  x  H  x  (0, T)  to  hi  x  x  (0,T).  Then  we  use  the 


explicit  representations  of  and  u,  equations  (6.19)  and  (6.28),  respectively,  to  show 
that  for  hxed  t  G  [0,T], 


lim 

£— »'0 


u^(x,  t)  —  u  f  j 


LP(0)3 


=  0. 


(6.36) 


Now  we  justify  Steps  Q-g. 

Proof  of  Step  ([^.  Let  V  G  U’(Vt)^,  where  1  <  p  <  cxo.  Then  we  hrst  show  that 

(Al)'^v(x)  =  xf(x)(Hf)"'n(a;)  for  all  n  G  N.  (6.37) 
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The  proof  is  by  induction  on  n.  The  formula  (6.37)  holds  for  n  =  1  by  the  dehnition  of  Af. 
Assume  that  it  holds  for  n  =  k.  Then  for  n  =  k  +  1, 


(AD^+Mx)  =  xl{x)Q  j  ^  ((A^)  V(x  +  (/,  0,  0))  -  (A^)^t;i(a;))  dl 

=  xl{x)Ci  j ^  {xl{x  +  (/,  0,  0))(Af)^'i;i(a:  +  (/,  0,  0)) 

-  xt{x){Af)’^vi{x))  dl.  (6.38) 

Note  that  since  x  lies  in  a  hber  if  and  only  if  a;  +  {I,  0,  0)  lies  in  the  same  hber,  then 
Xf{x  +  (Z,  0,  0))  =  Xf{x).  On  the  other  hand  (xf)^  =  xf,  thus  (6.38)  becomes 

=  xlix)  J  ((Af)  V(a:  +  (/,  0,  0))  -  (Af)  V(a:))  d/ 

=  Xf(a^)(Af)*^+^n(a:). 


Therefore  (6.37)  follows.  Since  (Af)"n  G  it  follows  from  Propositions  4.3  and  4.4  of 

Section  14.11  that 


xl{x){Ai)^v{x)  Xi{y){ATv{x). 


(6.39) 


Next  we  show  that 


(ATv{x)  ^  Xi{y){Airv{x). 


(6.40) 


To  see  this,  we  note  that  from  (6.17),  the  operator  (A^)*^  ,  n  G  N,  can  be  written  in  the 
following  form 


(AS)n  ^  ^  ^ 

where  the  operator  is  bounded  on  and  satishes 

ll^nll 

for  some  M  >  0  independent  of  £.  It  follows  that  for  hxed  n  G  N, 

MmeD^v  =  0,  in 
£— >0 


(6.41) 


(6.42) 


(6.43) 


and  thus  by  Proposition  |4.2[  the  sequence  {eDf^v)e>o  two-scale  converges  to  0.  Therefore  the 
result  follows  by  combining  (6.41),  (6.39),  and  (6.37). 


Now  we  recall  from  (|6.19)  that  u^{x,t)  is  given  by 


^2n+l 


M"(a:,  t)  =  u\x)  +  tv\x)  +  ^  (A")  V(a;)  +  ^  (A")  V(a:) 

^(2n)!  ^(2^  +  1)! 


(6.44) 
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Using  (6.40),  we  will  show  in  Section  6.2.1  that  for  G  /C, 

°°  i2n 


/*  j.'2n 

^  f2n 

Xf(2/)(^f)"M°(a;)  •  ij  {x,  y)  dydx,  (6.45) 

n=l  '  '' 


r  UL  +2n+l  ,  . 

lim  /  ^  (A^Yv^ix)  ■  ijj  (x,-)  dx 


^2n+l 

y]  77, — ^TT  Xi{y){AfTv^{x)  ■  Y  {x,  y)  dydx. 
n  JY  +  J-j! 


(6.46) 


It  follows  from  (6.45)  and  (6.46)  that  for  hxec 


u^{x,t)  u{x,y,t),  where  u  is  given  by  (6.28) 


t  G  [0, T],  as  £  — >•  0, 


Proof  of  Step  ([^.  We  can  see  from  (6.28)  that  u  G  C^([0,T];  L^(f2  x  U)^).  Then  by 
taking  the  second  time  derivative  of  both  sides  (6.28),  we  obtain 


^2n+l 


dMx, 2/, t)  =  Xi{y)  TTh  +  Xf{y)  Y1  (o  i  i^i 


=  xi{y)Mu° +  txYix) 

°°  y.2n 

+Xi{y)Aj2  7^  iyuYx)  +  Xf{y)AiJ2  (2n  +  l)\ 

n=l  ^  h  K  )■ 


°°  ^2n+l 


(6.47) 


From  (6.28)  and  the  dehnition  of  Af,  given  by  (6.30),  we  see  that 


°°  j.2n 


Aiu{x,yR)  -  7lf(M°  +  tn°)(a;)  =  Xf(?/)7lf^  — —  (7lf)V(a;) 

n=l  ^  ' 


^  y-2n+l 
n=l  ^  ' 


Thus  from  (6.47)  and  (6.48)  we  obtain  that 

dtu{x,y,t)  =  Aiu{x,y,t)  +  (xf(2/)  -  l)Ai{u°  +  tv^ix), 


(6.48) 


(6.49) 


and  hence  (6.29)  follows.  The  linear  operator  Af  is  bounded  on  LP{fl  x  U)^.  Thus  u  is  the 
unique  solut 


ion  of  (6.29) 


Proof  of  Step  (|^.  From  (|6.27[)  and  Proposition  |4.9[  we  obtain  that  for  hxed 
t  G  [0, T],  as  £  — >  0, 

^  j  u{-,y,t)dy  weakly  in  L^(f2)^. 
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By  definition  u^{x,t)  =  Jyu{x,y,t)dy,  thus  (6.32)  follows.  It  is  clear  that  (6.31)  follows 
from  integrating  (6.29)  over  Y. 


Proof  of  Step  (|^.  Dehne 

oo  5 
1=1  ^  ^ 


°°  f2n  °°  j-2n+l 

=  +  E  (2n  +  1)! 

n=l  '  ' 


(6.50) 


Combining  this  equation  with  (6.28)  gives  (6.33).  On  the  other  hand,  equation  (6.50)  implies 
that  w  G  C^([0,  T];  Thus  by  taking  the  second  time  derivative  of  both  sides  of  (6.50) 

gives 

°°  j.2n  j.2n+l 

(°°  j.2n  “  j.2n+l  \ 

^  M  ^  (&TT)!  j 

=  Af{u^  +  tv^){x)  +  Aiw{x,t).  (6.51) 


Note  that  from  (6.50)  it  is  easy  to  see  that  w(x,0)  =  0  and  dtw{x,0)  =  0.  Combining  this 
fact  with  (6.51),  equation  (6.34)  follows.  The  fact  that  Af  is  linear  and  bounded  on  Lp(0)^ 
implies  that  w  is  the  unique  solution  of  (6.34). 

Proof  of  Step  ([^.  Extend  xt  from  Y  to  by  periodicity.  Then  by  making  the 
substitution  y  =  ^  in  (6.28),  we  obtain 

OO  xStt. 

^(a:,  pt)  =  M°(a;) +fn°(a;) +  Xf(x)^^^  (klf)V(a;) 


+  xl{x)J2 


°°  ^2n+l 


(klf)V(x) 


^  (2n  +  l)! 


j-Zn 

=  u\x)  +  tv\x)  +  TTp 

n=l 

^  f2n+l 
n=l  '  ^ 


(6.52) 


where  in  the  last  equality  we  have  used  equation  (6.37). 
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Now  we  compute  the  difference  u^{x,t)  —u{x,  |,t)  using  equations  (6.19)  and  (6.52).  We 
see  that 


^  ^2n 


u%xR) -u(x,^R^  = 

n=l  ^  ^ 

^  j-2n+l 

((^'r-(An").V) 

71=1  ^  ' 


^  ^2n 


E 


S  (2n)! 


j-2n+l 

i^^n)u\x) +  J2  r2n  +  l)\ 

n=l  ' 


(6.53) 


where  in  the  last  equality  we  have  used  equation  (6.41).  By  taking  the  norm  in  (6.53) 


and  by  using  (6.42),  we  see  that 


u‘^{x,  t)  —  u  j 


°°  ^2n 


—  ^  (O  \\  ^'^\\'^^\\Lp{np 

LP{Q)3  ^  (2n)!  ^  ^ 

n=l  ^  ' 

^  j-2n+l 

+  ^.rniL.O 


^  (2n  +  l)! 


||lp(0)3 


71=1 


=  e  ^||m°||lp(o)3  coshA/Mt  +  ||n°||LP(n)3--^  sinh  a/M^^ 


thus  (6.36)  follows,  completing  the  proof. 


6.2.1  Proof  of  (6.45)  and  (6.46) 


In  this  section  we  prove  (6.45).  Equation  (6.46)  can  be  derived  similarly. 


We  begin  by  the  following  observation 


E/„  (ly  ^  i) 


71  =  1 


dx  <  oo. 


To  see  this,  we  use  Cauchy-Schwarz  inequality  to  obtain 
[  {ATu\x)  •  (x,  -)  dx  <  ||(kl^)VI 


i,2(n)3 


V'  (x,  ^ 


L2(n)3 


(6.54) 


(6.55) 


From  Part  (a)  of  Proposition  6.4,  the  operator  A^  is  uniformly  bounded  on  L^(f2)^.  Also,  it 
is  easy  to  see  that 

^(xA)  <  \mLHn-,c,MY)^)  ■=  (  [  sup  \^/J{x,y)\‘^dx\  . 

V  e/  L2(n)3  \.lny&Y  I 


We  use  these  two  facts  in  (6.55)  to  obtain 


(A^)V(x)  •  ip  (^x,  dx  <  ||w°|L2(f^)3  ||'0l|L2(O;CperT)3)>  (6-56) 
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for  some  M  >  0.  Therefore 

^2n 


n=l 


dx  <  ||m°I 


L2(Q)3 


L2{n-,Cper(Y)3) 


E 


^  (2n)! 


from  which  (|6.54)  follows. 

Now  from  (6.54)  and  by  using  Lebesgue’s  dominated  convergence  theorem,  it  is  straight¬ 


forward  to  show  that 

^  /-2n 


P  4-272  P  -f-OiTX 

/„  E  ^  ( V)”«“(x)  .  i,  (x,  f )  <i:.  =  g  ^  ( V)”«»(x)  .  i,  (x,  I) 


(6.57) 


For  n  G  N,  we  define 


SN,e  = 


N 


j.2n 


^  (2n)< 


72: 


•  'ip  dx. 


Then  using  (6.40)  we  see  that 

N 


x2n  r 

lim5Ar£  =  y^  /  ; 


Xf(l/)(^f)’"w  (a:)  ■i>{x,y)  dxdy, 


(6.58) 


and  hence 


lim  lim  SN,e  =  Y]  jttR  /  Xf(2/)(7ff)”'u°(a;)  •  -0  (x,  ?/)  dxd?/. 

^  (_2n)! 


(6.59) 


Below  we  will  show  that  the  order  of  the  limits  in  (6.59)  can  be  interchanged,  i.e.. 


lim  lim  SNe=  bm  limS'Afe. 

e^O  N^oo  ’  N^oo  e^O 


(6.60) 


Combining  this  with  (6.57)  we  obtain 

P  -4272 

l‘hy„E(^(^T«“W'V-(^.j)<i. 


^  ^2n 


E 


^  (2n)! 


Xi{y){ATu^{x)  ■  ^  {x,y)  dxdy. 


'OxY 


(6.61) 


Applying  arguments  similar  to  those  used  in  obtaining  (6.57),  we  can  show  that 

4272 


E 

n=l  ■ 


'OxY 


(2n)! 


Xi{y){ATu^{x)  ■  Ip  {x,y)  dxdy 


-  ^2n 


TTf  YyYiTu^{x)  ■  pj  (x,  y)  dxdy. 


'nxY  7^  (2^) 


(6.62) 
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From  (6.61)  and  (6.62),  the  result  (6.45)  follows. 

To  complete  the  proof,  it  remains  to  justify  (6.60).  It  is  sufficient  to  show  the  double 
sequence  is  Cauchy.  So  assume  that  iV,  L  G  N  such  that  N  >  L.  Then 


< 


< 


n=L+l 

N 

E 

n=L+l 


{2n)\ 

e 

(2n)! 


in 


■  4)  (^x,  -  j  dx 
{A^Yu^{x)  ■  4j  ^x, 


N 


h°llL2(o)3 


dx 


(6.63) 


n=L+l 
N 


where  (6.56)  was  used  in  the  last  step.  We  note  that  the  term  X]n=L+i  i'^  (6.63)  can 

be  made  arbitrarily  small  by  choosing  large  values  of  N  and  L.  We  conclude  that  for  given 
C  >  0,  there  exists  a  positive  integer  K{Q  such  that  for  N,L  >  K{()  and  all  £  >  0, 


\SN,e  -  SL,e\  <  C 


(6.64) 


From  (|6.58)  and  (6.64),  and  by  using  Lemma  6.6  below,  it  follows  that  the  double  sequence 
{SN,e)  is  Cauchy. 

Lemma  6.6.  Let  {an,k)  be  a  double  sequence  in  d  eN,  such  that 

(a)  For  each  n  eN, 

lim  a^i^k 

k^oo 

(b)  Given  C  >  0,  there  exists  a  positive  integer  N  =  N{()  such  that  for  n,l  >  N  and  all 
kEN, 


\an,k  -  ai^k\  <  C-  (6.65) 

Then  the  double  sequence  {an,k)  is  Cauchy,  and  hence  convergent. 

Proof.  Let  C  >  0  ^ind  assume  that  N  E  N  satisfies  Part  (b).  Then  consider  the  sequence 
{ciN,k)k£N-  It  follows  from  Part  (a)  that  this  sequence  is  convergent,  and  hence  Cauchy.  Thus 
there  exists  a  positive  integer  K  =  K{N,  ()  such  that  for  k,m  >  K, 


\(iN,k  —  0'N,m\  <  C- 


Let  J  =  max{iV,  iP}.  Then  from  (6.65)  and  (6.66)  we  obtain  that  for 
n,  l,k,m  >  J, 


(6.66) 


—  ^N^k\  \^N,k  ^N,m\  \^N,m 

<  3C, 


and  therefore  the  double  sequence  (a^^fc)  is  Cauchy. 


□ 
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